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ABSTRACT 



Context. The evolution of intermediate and low-mass stars on the asymptotic giant branch is dominated by their strong dust-driven 
winds. More massive stars evolve into red supergiants with a similar envelope structure and strong wind. These stellar winds are a 
prime source for the chemical enrichment of the interstellar medium. 

Aims. We aim to (1) set up simple and general analytical expressions to estimate mass-loss rates of evolved stars, and (2) from those 
(~| [ calculate estimates for the mass-loss rates of the asymptotic giant branch, red supergiant, and yellow hypergiant stars in our galactic 

sample. 

I ■ Methods. The rotationally excited lines of carbon monoxide (CO) are a classic and very robust diagnostic in the study of circumstellar 

O ' envelopes. When sampling different layers of the circumstellar envelope, observations of these molecular lines lead to detailed pro- 

files of kinetic temperature, expansion velocity, and density. A state-of-the-art, nonlocal thermal equilibrium, and co-moving frame 
C/3 , radiative transfer code that predicts CO line intensities in the circumstellar envelopes of late-type stars is used in deriving relations 

between stellar and molecular-line parameters, on the one hand, and mass-loss rate, on the other These expressions are applied to our 
extensive CO data set to estimate the mass-loss rates of 47 sample stars. 

Results. We present analytical expressions for estimating the mass-loss rates of evolved stellar objects for 8 rotational transitions of 
the CO molecule and thencompare our results to those of previous studies. Our expressions account for line saturation and resolving 
of the envelope, thereby allowing accurate determination of very high mass-loss rates. We argue that, for estimates based on a single 
rotational line, the C0(2-l) transition provides the most reliable mass-loss rate. The mass-loss rates calculated for the asympotic giant 
branch stars range from 4 x 10"^ Mq yr"' up to 8 x 10"^ Mq yr"'. For red supergiants they reach values between 2 x 10"' Mq yr"' 
and 3 x 10"* Mq yr"' . The estimates for the set of CO transitions allow time variability to be identified in the mass-loss rate. Possible 
mass-loss-rate variability is traced for 7 of the sample stars. We find a clear relation between the pulsation periods of the asympotic 
00 ' giant branch stars and their derived mass-loss rates, with a levelling off at ~3 x 10"^ Mq yr"' for periods exceeding 850 days. 

' Conclusions. 



Key words, stars: AGB and Post-AGB - stars: Supergiants - stars: mass loss 



1. Introduction stars, they are typified by circumstellar envelopes of dusty and 

^ molecular material created by strong winds. The mechanisms 

■ Stars of low and mtermediate masses (1 M^<M^< 9M^) expel ^^^^^^ den se outflows of b oth types of objects are likely 
a large part of their outer layers m a dust-driven wind during the ^^^^^^^ (I.Tosselin & PlezL ^07), since the central stai-s 

later stages of their evolution, when ascending the asymptotic ^^^^ obtained other physical and chemical properties after go- 
giant branch (AGB). The expulsion of stellar material creates ^^^^^^^ somewhat different phases of evolution, 
a cool circumstellar envelope (CSE) containing dust grains and 

molecular gas phase species. The rates at which mass is lost (M) We want to investigate the physical and chemical proper- 
vary from ~10"^ Mq yr"' up to ~10"'* Mq yr"' . A violent ejec- ties of these CSEs by analysing molecular-line observations. The 
tion of the circumstellar envelope ends the AGB phase and ini- most important species for such investigations is carbon monox- 
tiates further evolution into a protoplanetary nebula (PPN) and ide, CO. More specifically, observations of lines resulting from 
later int o a planetary nebula (PN) when the ejected material is transitions from states with angular momentum quantum num- 
ionised dHabing & Olofssonll2003l) . ber 7 to 7 - 1 can be used to derive density, velocity and kinetic 
High-mass stars (M* > 9 Mq) do not ascend the AGB, temperature in the envelope. Previous surveys have focussed on 
but evolve into supergiants and hypergiants. Similar to AGB the analysis of only a few low-excitation rotational transitions 
of CO (7 = 1-0, 2-1, 3 - 2), or have used simple ana- 
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eters of these tratisiti otis (iKnapp et all I1982L iKnapp & MorrM 
1985; 'Olofsson et al.', '1 993l:lLoup et al.Lll993l:lNeriet all 1199a 
Ramstedt et al., 2008). Since we have access to an extensive data 
set covering both low and high-/ molecular lines, we want to 
provide mass-loss estimators based on multiple CO lines, includ- 
ing those of high excitation levels. These analytical expressions 
will link stellar and CO-line parameters to the (variable) mass- 
loss rates of the central stars. 

The presence of both '^CO and '^CO line transitions in the 
data set allows us to estimate the '^C/'-'C isotope ratio for many 
stars in the sample. This abundance ratio provides a measure for 
the chemical evolution and hence of the evolutionary state of 
the stars. In this paper we give first order estimates for '^C/'^C, 
derived from line intensity ratios. 

The sample of evolved stars and the CO-line observations 
are presented in Sect.|2]and in the (online) appendix. We discuss 
some of the sample stars individually based on the observations. 
Sect. [3] deals with the radiative transfer analysis and the param- 
eter study that will lead to the construction of analytical expres- 
sions to estimate mass-loss rates. The constructed formalism is 
then compared to those presented in the literature. In Sect.|4]and 
Appendix we discuss the determination of the basic stellar 
parameters needed to derive mass-loss rates. The results of our 
study in terms of M, '^C/'^^C and wind driving efficiency are 
presented in Sect. |5] Conclusions are given in Sect. |6] 

2. Observations 

2.1. The sample 

The data for our sample of 69 galactic objects have been assem- 
bled over many years and mainly consists of targets originally 
selected as potential candidates to be included in guaranteed 
time key programs (GTKP) of the Herschel-HIFI and Herschel- 
PACS instruments. Different evolutionary and chemical types 
are covered, but there is a rather strong bias towards oxygen- 
rich (C/0< 1) AGB stars. Several categories of AGB stars, dif- 
fering in pulsational and mass-loss properties, are represented 
in the sample, i.e. Mira-variables, semi-regular (SR) variables 
and OH/IR-type stars. SRs and Miras have low to intermediate 
mass-loss rates (10^*^ - 1O"''M0 yr"'), while OH/IR stars have 
very dusty and optically opaque envelopes, formed by intense 
mass loss that can reach up to a few 10""^ Mq yr"'. Typical for 
these objects are the strong OH masering at 1612 MHz and a 
very large infrared excess. Miras have very regular pulsations, 
with a fixed period of a few 100 days (average ~400 days) 
and large amplitudes (>2.5mag in V-band). Semi-regulars, on 
the other hand, have smaller amplitudes and can exhibit ureg- 
ularities in their pulsations. Their pulsation per iods can range 
from 20 up to 2000 days (GCVS 'Samus_eLaLl l200l . OH/IR 
stars have pulsation periods ranging from a few 100 days up 
to more than 1000 days with an average of ~1000 days and 
they could be considered th e evolutionary successors of Miras 
(IVassiliadis & Woodl Il993h . Their thick envelopes, obscuring 
the central star in the optical, are produced by the so-called 
super wind phases (Iben & Renzini, 1983; Vassiliadis & Wood, 
Il993h . which occur during the late parts of the quiescent hy- 
drogen burning phase, just before the central star goes through 
a thermal pulse (TP). During these phases, the pulsation period 
increases with increasing mass loss (up to a few 10"^ M© yr"') 
and luminosity. The enhanced mass loss leads to a very optically 
thick and dusty CSE, characterising an OH/IR star When a su- 
perwind phase is halted, the pulsation period decreases again. 



the expelled material diff uses outwards, and a less in tense mass- 
loss process is initiated (IVassiliadis & Woodi Il993h . The com- 
bination of these factors can lead to the central star being again 
visible in the optical, and being again observed and classified as 
a Mira or SR variable. This scenario is more plausible for the 
less massive stars (M* < 2.5 Mq), as the higher-mass stars (M* 
> 2.5 Mq) experience relatively modest variations in their high 
mass-loss rates, lumino sities, and pulsation period s due to the 
TPs (see Figs. 3 to 9 in IVassihadis & Woodl 119931) . Low-mass 
stars experience only very short periods of enhanced mass loss. 
Moreover, during these periods the absolute value of the mass- 
loss rate is significantly lower than for their massive counter- 
parts. Therefore, the chance to observe them in their superwind 
phases, i.e. as OH/IR, is far smaller than for the higher-mass 
stars. 

Supergiants and hypergiants and some post- AGB objects are 
also included in the sample. The latter fit in this study as the 
progeny of AGB stars. Studying their extended envelopes will 
lead to an improved knowledge of the envelopes of AGB -type 
stars and the late stages of the AGB evolution. Two young stellar 
objects (YSOs) — AFGL5502 and the Gomez Nebula — were 
observed together with the sample of evolved stars. The obtained 
data on these YSOs were not yet published and are presented in 
this paper. We will not further discuss these objects, since the 
focus of this paper is on evolved stars. 

2.2. The observations 

The data set presented in this paper consists of observations of 
multiple rotationally excited lines of both '^CO and '^CO in 69 
star s. The bu l k of t he presented lines were first published by 
Kem per et al.l (l2003h . Transitions from J — 1-0 up to 7 = 7- 6 
for '^CO and from 7 = 2- 1 up to 7 = 6-5 for '^CO are covered. 
Since the high-/ lines have formation regions deeper within the 
CSEs, this multitude of rotational lines can sample a much larger 
part of the CSE than only the low-excitation transitions. Also, 
because of slightly different molecular properties, the '^CO-data 
probe regions somewhat different from those sampled with the 
'^CO -rotational lines. Moreover, since the abundance of '^CO 
is lower than that of '^CO, the circumstellar layers are opti- 
cally thin(ner) for the rotational lines of the former, implying 
that CO -rotational lines provide good diagnostics for the enve- 
lope's density structure. In this respect, we should be able to ob- 
tain a more complete view on the structure and properties of the 
envelopes and the mass-loss history of the objects, than was pos- 
sible before with single-dish data only sampling low-excitation 
i^CO-lines. 

All previously unpublished data in our sample were obtained 
with APEXQ (Atacama Pathfinder Experiment) and JCM10 
(James Clerk Maxwell Telescope). Other data were obtained via 
private communication or retrieved from the literature. 

APEX is a 12m single-dish telescope with a frequency range 
from 210 up to 1500GHz located at Llano Chajnantor, Chile. 
CO-line data were obtained with three heterodyne SIS-receivers 
mounted on the Nasmyth-A focus: APEX-2A, FLASH-I and 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 

- The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 
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Table 1. Frequency ranges, half power beam widths Ob, main 
beam efficiencies 77MB, and observable CO rotational excitation 
lines for APEX and JCMT instruments. 



Instrument 


Frequency 




MB 


Observable 






(GHz) 


(arc sec) 




CO lines 




APEX 












APEX-2A 


279 _ 381 


17.3 


0.73 


'^CO, ''CO (3 - 


2) 


FLASH-I 


430 - 492 


13.3 


0.60 


'^CO, "CO (4- 


3) 


FLASH-II 


780 - 887 


7.7 


0.43 


'^CO (7 - 6) 




JCMT 












A 


211-279 


20 


0.69 


'-CO, '''CO (2 - 


1) 


B 


315-375 


14 


0.63 


'-CO,'3CO (3- 


2) 


HARP 


325 - 375 


14 


0.63 


'-C0, ''CO (3 - 


2) 


RxW(C) 


430-510 


11 


0.52 


'-CO, '''CO (4 - 


3) 


RxW(D) 


626 - 710 


8 


0.30 


'^CO, ''CO (6 - 


5) 


E 


790 - 840 


7 


0.25 


'^CO (7 - 6) 





FLASH-II. Instrument specifics are listed in Table [U together 
with the CO-lines observable in the respective frequency ranges. 
All observations were performed in beam-switching mode. 

Reduction of the APEX data was done with CLASS, part of 
the GILDAS-package. The process mainly consisted of combin- 
ing different scans of the same molecular line towards one object 
into a single spectrum and subsequently removing baselines and 
spikes in the obtained spectra. Correcting the data for telescope 
efficiencies was carried out by the APEX online calibrator. 

JCMT data were obtained with heterodyne receivers A, 
B, HARP, RxW (bands C and D), and E in beam-switching 
mode. Data reduction was performed with the reduction pack- 
age SPECX. Spectra of standard stars taken during the observ- 
ing runs were compared to the standard spectra available on the 
web site of JCMT. In case of large deviations in the measured 
line intensities of these standard stars (> 10%), the scientific- 
programme data obtained around the time of the measurement 
of these standard stars were corrected for these discrepancies. 
In case of multiple available standard spectra with other mag- 
nitudes of deviation, the correction factor was determined via 
linear interpolation in time. As for APEX data, most corrections 
and conversions of the data are performed online. The correction 
for the main-beam efficiency t/mb had to be done manually in the 
course of the reduction process using the values given in Table[T] 
;7mb allows transforming antenna temperatures, Ta, into main- 
beam-brightness temperatures, Tmb- The latter is the equivalent 
of the brightness temperature of measurements performed with 
a perfect antenna outside the earth's atmosphere. The main ad- 
vantage of a Tmb -based temperature scale is that the data are 
no longer dependent of any of the instrumental properties, apart 
from the beam width. All data in this paper are presented on the 
JiviB-scale. 

The uncertainties on the JCMT data, so-called absolute er- 
rors, are fairly well established for all receivers (Kemper et al., 
120031) . For APEX this is not so clear. No standard spectra were 
available to check the performance of the instruments at the time 
of the execution of the scientific programme and no standard val- 
ues for the uncertainties on the data were found in the literature. 
[Ramstedt et al. (2006) mention an uncertainty of 20% on the 
absolute intensity scale for the 7 = 3-2 transition, but give no 
reference for this number 
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Fig. 1. '^CO and '^CO data of IKTau, obtained with JCMT and 
APEX. It is clear from the graphs that the '^CO lines are more 
parabolic in shape than the '^CO lines. This effect is caused by 
the lower abundance of '^CO, leading to optically thinner line 
profiles. 



2.3. Resu/fs 

Tables lA.ll and IA.3I in Appendix |A] (available online) contain 
main-beam-brightness temperatures at the centre of the line pro- 
file (Tmb.c), velocity integrated main-beam intensities (/mb), and 
expansion velocities (voo) for all CO-data of the sample stars. 
Among the high-quality CO-data for over fifty stars, there are 
'^CO data for 29 targets. 

Some of the observed lines towards AGB or RSG stars 
in the sample were not detected, e.g. '^CO (3 - 2) for 
ffSco (=IRAS 16262-2619), and some others were too noisy 
or heavily contaminated by interstellar lines, e.g. V1360Aql 
(=IRAS 18432-0149). In these cases, the data are presented in 
the appendix, but are not used in the further analysis of the sam- 
ple. The respective targets are listed in Table |5] together with 
Post-AGB objects and YSOs. 

Characterisations of the sample targets in terms of chemical 
type (oxygen-rich, carbon-rich or S-type) and pulsational and/or 
evolutionary type (e.g. Mira, SR, OH/IR) are listed in Table|6] 

2.3.1. Line profile diagnostics 

Because of its high stability against photodissociation and its 
large molecular abundance, carbon monoxide exists throughout 
the largest part of the CSE. Rotational transitions of CO are eas- 
ily excited in the cool envelopes and can therefore trace many 
properties of all different layers, e.g. density and temperature. 
These properties are linked to the optical thickness of the enve- 
lope layers and the mass-loss rate of the central star. 

In case of optically thick, unresolved envelopes, the ob- 
served line profiles have parabolic shapes. Changing to the 
optically thin and/or spatially resolved case, the molecular 
lines are more flat-topp ed or even display two-horned profiles 
(iKnapp & Morrisl|l985l) . Considering this simple rule, the shape 
of the line profiles is a useful first diagnostic for the (density) 
properties and the geometric extent of a star's CSE relative to the 
telescope beam. The data of, e.g., IKTau shown in Fig.[T]clearly 
reflect that the envelope layers are optically thinner for '^CO 
than for '^CO, since the line profiles of the latter are consis- 
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Fig. 2. Panels a to c show the multiple superimposed peaks that are easily detected in the CO emission lines of R Cas, V Hya, 
and VY CMa. Panels d and e show the two separate velocity components seen in the CO lines of U Ant and S Set. These are clear 
evidence for detached shells ai^ound these targets. It has also been shown that R Scl (panel f) has a detached shell. See text for further 
comments. 



tently more parabolic in shap e. The mass loss of CW Leo is very 
strong (~2 x 10 ^''Mq vr~': iMauron & Hugginsl l2000l) . caus- 
ing optically thick CO emission lines. These are, however, flat- 
tened due to the larg e angular size of the envelope (~200 arcsec; 
iMauron & Hugginsl I200Q) . causing resolution by the telescope 



beams with half power beam widths between 7 and 20 arcsec 
(see Table[T]l. 

A first analysis of the data co nsists in fitting a so-called soft 
parabola dOlofsson et all Il993h to the observed line profiles. 
Every rotational line is specified by several line parameters: the 
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Fig. 3. The CO line profiles measured towards the oxygen-rich 
SRb variable EP Aqr bear clear evidence of a composite circum- 
stellar environment. See Sect. l2.3.2l for further details. 



main-beam-brightness temperature at the line centre, Tmb,c, the 
velocity at the line centre, i.e. the velocity of the star with re- 
spect to the Local Standard of Rest, vlsr, and the half width of 
the line profile, i.e. the expansion velocity of the CSE in the out- 
ermost regions where the studied molecule is present, Voo. These 
parameters are obtained by fitting the data with the soft parabola 
line profile function, given in Eq.[T] where yS is a measure for the 
shape of the profile function: 



T(v) 



' MB,c 



1 



Vlsr 



13/2 



■ Vlsr I < Vc 



(1) 



0, 



■ Vlsr I > 



When p - 2, the line profile has a parabolic shape, represen- 
tative for optically thick lines, observed towards spatially un- 
resolved CSEs. Smaller, positive values for p lead to more flat 
topped profiles. Negative fi values can be used to fit two-horned 
profiles observed towards optically thin envelopes. The best fit to 
the line profile was determined through minimising the total ab- 
solute difference between the data and the soft-parabola fit with 
the IDL-routine AMOEBA. 

Figs. lA.l] and lA.2l show all '^CO and '^CO observations and 
the soft-parabola fits to these data. The y6-values per line profile 
are listed in Tables lA.ll and IA.3I The parameter fi provides in- 
formation on the optical thickness of the envelope (and hence 
on M) and the resolving power of the instrument. Its value will 
be used to improve the accuracy of the mass-loss determination 
(see Sect. [Ml- 

2.3.2. Discussion on individual targets 

If strong deviations from the soft-parabola fit are found, this fit- 
ting procedure could reveal significant deviations from the as- 
sumptions of, e.g., spherical symmetry or constant M, since de- 
tached shells, variability in the mass loss, vast clumps in the en- 
velope, or jets can affect the symmetry and the overall shape of 
the CO lines. Some of the targets in the sample and their devia- 
tions from the soft-parabola fits are discussed individually in this 
section. 

It is obvious from Figs. lA.lj and lA. 21 that indeed not all pro- 
files in the data set can be fitted with the simple type of line 



profile function presented in Eq. [T] R Cas (Mira) and V Hya 
(SRa) have line profiles exhibiting superimposed peaks — see 
Fig.[2 | — that could be link e d to bipolar out flows ( Olofsson et all 



ll993riSahai etali l2003h . lOlofsson et alJ (fT99 3) suggest that 
these bipolar outflows could be related to th e possible presence 
of a binary companion in case of both targets. [Sahai et al] (l2003h 
suggest that V Hya could already be transitioning to the Post- 
AGB phase. The high- velocity bipolar outflow would be linked 

to this transition process. 

U Ant and S Set are discussed by lOlofsson et alJ (Il993h as 
targets with convincing evidence in the CO data for highly 
episodic mass loss, possibly due to a thermal pulse. The shapes 
of their line profiles — see Fig. |2] — indeed imply the presence 
of a detached shell. The inner parts of the outflow have terminal 
velocities of respectively 7kms"' and 8kms"', while the out- 
flow velocities of the outer parts of the envelopes reach 1 6 km s ' 
and 20 km s ^ The carbon-r i ch se mi-regular R Scl also has a de- 
tached shell (lOlofsson et all [19961) . albeit not directly visible in 
the CO line profiles. 

The CO line profiles of the SRb vaiiable EPAqr (Fig. [3]) 
reveal the composite nature of its circumstellar environments. 
Except for the low-S/N observation of '^^CO(4-3), all mea- 
sured line profiles towards EPAqr exhibit a composite profile 
with (a) a broad, Iow-Tmb profile and (b) a very narrow, high- 
T^MB profile centrally superposed on this broad profile. Only 
very little is known about the origin of these types of line pro- 
files. The broad plateau emission has been suggested to origi- 
nate from episodic mass loss, bipolar out flows, or circumstel - 
lar disks (KnapD et al., 1998; Kerschbau m & OlofssonL Il999l) . 
[Winters et al.l (1200 7) mention that no obvious departure from 
circular symmetry can be seen, but that there is evidence for a 
ringed structure in the '^CO(2-l) map implying variation of M 
in time. 

RHy a is well-known for its declin ing per iod and mass- 
loss rate (Dec in et"aLL l2008Natur eexlabai) . which [Wood & Zarol 
(.1981.) attributed to a possible recent thermal pulse. The de- 
tached shell that has been dete cted at 60 nm can be explained 
by a slowing down of the wind. lDecin et alJ (l2008Natureexlab3) 
have modelled the wind of RHya in detail, checking mod- 
els produced with the radiative transfer code GASTRoNOoM 
against data of both rotational and vibration-rotation transitions 
of CO. WXPsc is a second object for which M could be con- 
sidered to be under influence of TPs. This extreme oxygen- 
rich star cou ld currently be going through a superwind phase 
(iDecin et al 1^007), causing the large infrared excess. For both 
targets, however, these presumed deviations from constancy in 
mass-loss rate are not directly visible in the rotational CO line 
profiles. 

The triply peaked low-7 transition profiles in both '^CO and 
'^CO towards the supergiant VYCMa — see Fig.[2| — are most 
probably a superposition of an optically thin component (red- 
wing and blue-wing peaks) and an op tically thick one (middle 
peak), as shown bv lZiurvs et al.l (l2007i) . For the higher-7 transi- 
tions (7 = 4-3 and up), the extra parabolic feature is no longer 
clearly present. This, in combination with the asymmetry of the 
profile, led to the assumption of a variable mass-loss history as 
described by Decin et al. (2006). NMLCyg and aOri both have 
spiked line profiles, somewhat similar to those of V Y CMa. In 
the case of NML Cyg, the lines could be contaminated with inter- 
stellar CO detection. The line profiles of IRC -1-10420 and yU Cep 
exhibit asymmetries, with a blue-wing feature superimposed on 
a more or less parabolic profile. The extra feature is far more 
pronounced for // Cep and is also further out to the blue part of 
the line profile. 
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For some of the more evolved objects in our sample, i.e. 
Post-AGB stars, we also see deviations from the soft-parabola 
shape. The ^^CO(3 - 2) line profile towards GSC 08608-00509 
exhibits a sharp red-wing peak. For the Cotton Candy Nebula, 
we see a similar sharp feature in the blue wing. Since we have 
only few observations for these targets, it is hard to draw well- 
founded conclusions about the properties or geometry of their 
extended envelopes. 

All three line profiles towards the Gomez Nebula are 
two-horned, supporting the hypothesis of a CO emitting disk 
in keplerian rotation, as mentioned by iBuiarrabal et al.l (l2008l 
their Fig. 3). The CO line profiles towards the other young 
stellar object, AFGL5502, bear evidence of asymmetry, since 
the higher-excitation lines have superimposed features in their 
blue wings. 

We note here that TCet was not included in the sample of 
red s upergiants, contrary to what is often found in literature 
(e.g. ISloan&Pric5ll998l) . but in the sample of AGB stars. The 
low luminosity and the dust features in the ISO spectru m make 
classification of T Cet as AG B star more plausible CSpeck et al.L 
l2000HVerhoelst et al.Ll2009l) . 

3. Radiative transfer analysis 

Physical information on the envelope structure of evolved stars 
can be extracted from the observed line fluxes using a non-local 
thermodynamic equilibrium (non-LTE) radiative transfer code 
(see Sect. l3.1l for a short description). The availability of a large 
'^CO and '^CO line flux data base, as presented in Sect. 12.31 
off'ers interesting possibilities for statistical analyses. However, 
this large database also forces us to construct simplified ana- 
lytical expressions relating the mass-loss rate and the observed 
integrated CO-line intensities to deduce the mass-loss history 
of each target. We therefore derived new mass-loss rate formu- 
lae based on a grid computation. Similar formalisms that have 
been presented in the literature are summarised in Sect. 13.21 The 
grid construction and the resulting analytical expressions are dis- 
cussed in Sect. 13.31 and 13.41 A discussion on the quality of the 
formalism can be found in Sect. 



3.1. Radiative transfer model 

The structure of circumstellar envelopes is thought to be quite 
complex, due to the possible presence of shocks, bipolar out- 
flows, or episodic mass-loss events. The complex heating and 
cooling processes due to excitation of molecules, and the un- 
known grain-size distributions contribute to the complexity of 
the envelopes. The often taken approach to model these en- 
velope structures is therefore to choose a simplified approach 
with reasonable assumptions. We use the non-LTE radiative 
code GA STRoNOoM, w hich was presented and discussed by 
iDecin et al.1 (1200a l2007l l2008Nafiireexlabal) an d benchmarked 
accord ing to the method described by Ivan Zadelhoff et al.i 
(I2002h . considering a simple 2-level molecule (their cases la 
and lb). The CSE is assumed to be spherically symmetric and 
formed by a constant mass-loss rate. The code calculates the 
kinetic temperature and velocity structure in the shell by solv- 
ing the equations of motion of gas and dust and the energy bal- 
ance simultaneously lJusttanont et al.l (|1994|) . Adopting a value 
for the gas mass-loss rate, the dust-to-gas ratio is adjusted 
until the required terminal velocity is obtained (Decinetal., 
|2006|) . The statistical equilibrium equations and the radiative 
transfer equation are solved in the co-moving frame using 



the Approximated Newton- Raphson operator as developed by 
ISchonberg & Hemp e'('l986'). Finafly, a ray tracmg program uses 
the computed non-LTE level populations to calculate the line 
profiles . For a full d escription of the GASTRoNOoM-code, we 
refer to lDecin et all (12006^) . 

As describe d by iDecin et all (l2006l l2007l 
2008Natureexlaba), a detailed radiative transfer modelling 
of one target often results in the computation of many (> 10^) 
models. Some prior knowledge on the input parameters helps 
to restrict the extensive parameter space. The main input 
parameters for this research are the stellar temperature T*, 
the stellar radius the distance d, the CO fractional abun- 
dance /h,co = CO/H, the '^C/'^'C isotope abundance ratio, 
the terminal velocity Voo, the mass-loss rate M, and the dust 
condensation radius /dinner- When each of these eight parameters 
would only be scanned over five diff'erent values, almost 400 000 
models should be computed. Performing an in-depth line profile 
analysis for each target presented in this study is therefore 
too time-consuming and beyond the scope of this paper. The 
mass-loss rates of some twenty targets of the current sample 
will be individually modelled in a forthcoming paper (Lombaert 
et al., in prep.). 

Before describing the grid computation in Sect. 13.31 we give 
a short overview of mass-loss rate formulae presented in the lit- 
erature. This overview is for sure not complete, but provides a 
few often applied formulae, relevant to this study. These for- 
malisms are also compared with the one we present in this paper 

3.2. Comparison witli otiier studies 

Several theoretical mass-loss rate formulae based on observed 
line intensities or fluxes wer e already presented in the litera- 
ture. iKnapp & Morris I (Em derived a formula to estimate M 
in terms of the CO(7 = 1 - Oj3 line intensity. Their formula 
relates the mass-loss rate to an assumed CO abundance and a 
number of easily determined observables: the main beam tem- 
perature Tmb.c at the line centre (in K), the terminal velocity Voo 
(in kms"'), and distance d (in pc). A distinction was made be- 
tween optically thick and optically thin envelopes. The model 
for the CO emission assumed a spherically symmetric enve- 
lope expanding at constant velocity. The kinetic temperature 
profile through out the envelope was a ssumed in all cases to be 
that derived by iKwan & Linkd (11982) for CWLeo. The statis- 
tical equilibr ium equations were so lved for the lowest 15 rota- 
tional levels. Olofsson et aTl ([T99I transformed the formula of 
IKnapp & Morrisi (.19851) for an optically thick envelope to the 
slightly more general form 



M - S7v 10-20 

Mkm - 5.7 X 10 0.85 



Mo/yr, 



(2) 



where 0b is the beamwidth of the instrument (in arcsec) used for 
the observation of the line under study. The factor s{J) depends 
on the rotational transition J ^ J - I. For the / = 1 - tran- 
sition s(J) - 1, an d for the 7-2-1 transition s(J) - 0.5 (see 
the discussion in Olofs son et al.L ll993). In this equation — valid 
in the optically thick case — the CO envel ope size was fixed at 
3 X 10 ^^ cm, appropriate for high-M targets. fSchoier & OlofssonI 
('2001') compared mass-loss rates derived using their Monte 
Carlo non-LTE radiati ve transfer model to fl iose estimated us- 
ing the formulae of IKnapp & MorriH (Il985h . They found that 
the formula underestimated the mass-loss rates when compared 

^ When no isotope number is mentioned, CO is equivalent to '^CO. 
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CO(I-O) C0(2-1) C0(3-2) C0(4-3) 




-15 -10 -5 5 10 15 -15 -10 -5 5 10 15 -15 -10 -5 5 10 15 -15 -10 -5 5 10 15 

Velocity [km/s] Velocity [km/s] Velocity [km/s] Velocity [km/s] 



Fig.4. Line profile predictions for =2000K, ^SOOOLq, t/ = 300pc, /h,co = 1 x 10 ^ 12q13(;^^ jq^ /?i„nei = 3/?*, M= 1 x 
10^' Mq yr Vco = lOkms"' for abeam width of 10". The soft-parabola fit is plotted in dashed grey /light blue lines. For the high-/ 
excitation lines one can see both the stellar and dust continuum and the acceleration zone (the narrow peaks around the line centres). 



to those obtained from their radiative transfer analysis on aver- 
age by a factor of about four, and the discrepancies were found 
to inc rease for lower mass-loss rate objects ( Schoier & Olofsson, 
120011 see their Fig. 8). In the case of optically thin envelopes, the 
exact envelope size matters in estimating M, so Eq.|2l is indeed of 
lesser value. A correction for the envelope size (see iNeri et afl 
[1998 1 and references therein) would decrease the discrepancy 
found by Schoier & Olofsson (2001). However, considering the 
si mphcity of Eq.lH th e agreement is still remarkably good. 

iLoup et all (Il993h calculated the envelope size Rco and the 
mass-loss rate Ml by numerically solving the system of Eqs. [3] 
andg] 



Ml 
Ml 



2 X 103(2/h,co)»-85 

\7.3xl0if'/ \10/ 



xF-\Rco) Moyr-i 



1/1.45 



(3) 
(4) 



4x 10- 
2/h,co 



_4x 1/1.16 



Moyr 



They assumed /h,co of 2.5 x lO"'' for oxygen-rich stars and 
5 X 10 for both carbon-rich and S-type stars. The function 
F{Rco) introduces a correction factor for the envelope size and 
is normalised to unity for Rco - 3 x 10'^ cm. 



iRamstedt et alj (l2008l) derived a formula relating the mass- 
loss rate to the integrated intensity /co,y of the '^CO (7 
7-1) lines, the terminal velocity Voo, the abundance /h,co, 
the distance d to the object, and the telescope's beam size 6\,. 
Using the Monte Carlo r adiative transfer code described by 
ISchoier & OlofssonI (1200 ll) . the estimator was derived from a 
grid of 60 models varying in M (1, 3, 10, 30, and IOOxIO^^Mq 
yr-i), Voo (5, 10, 15, and 20kms-i), and /h,co (1, 3, and 
lOxlO"'*). The distance was taken to be lOOOpc. For each 
model, the velocity integrated intensities in CO (7 = 1 - 0, 2 - 
1, 3-2, 4-3) were calculated for a 20 m telescope (with cor- 
responding beam sizes of 33", 16.5", 11", and 8.5"). They use 
the s o-called /i-parameter to des cribe the free parameters of the 
dust (ISchoier & OlofssonI 1200 ih : 



h = 



/ 2.0gcm-n / 0.05 i^m 



(5) 



where *P is the dust-to-gas mass ratio, and pd and are the av- 
erage density and the radius of an individual dust grain, respec- 
tively. The value of h was assumed to be 0.2 for M in the range 
up to 3 X 10-^ Ma yr-i, 0.5 for 1 -3 x lO-'-Mo yr"', and 1.5 for 
10"^ Mq yr"'. The stellar temperature T*, stellar radius R^, or 
inner dust condensation radius Thinner were not specified. An esti- 
mate for the mass-loss rate was determined from a minimisation 



E. De Beck et al.: CO line survey on evolved stars and M-formuIae 



Table 2. Values for the input parameters used in the grid compu- 
tations. 



Param. 


Grid values 


Param. 


Grid values 




2000, 2500, 3000 K 


dinner 


3, 10, 30 iJ* 




3000, 10000, 3OOOOL0 


'^C/'^C 


10, 30, 100 


d 


300, 1000, 3000 pc 


M 


1, 3, 10, 30, 


Voo 


10, 15, 20kms-' 




100, 300, 1 000 


/h,co 


1, 3, 5x10-'' 




xlO-^Mo yr-' 


e 


10, 20, 30" 







procedure applied to Eq. |6l 



Mr = iy(/co,y0b«?VlO'')"^vlYH,coM0yr 



(6) 



with /co,y in K km s 



in km s , c/ in pc, and Oh in arcsec. A 



comparison between their input mass-loss rates and the values 
found from Eq. |6] was not shown, nor were estimates of the 
uncertainties on the derived M-values provided. The comparison 
between mass-loss rates derived using Eq.|6]and mass-loss rates 
derived using their Monte Carlo code for a sample of C-, Al- 
and S-stars shows quite a good agreement. Close inspection, 
however, reveals that for low-/ transitions (J - 1 - 0, 2 - 1) 
Eq.|6]tends to underestimate the mass-loss rate for M < IQ^^ Mq 
yr"' and to overestimate it for M > IO^^Mq yr"' (see Fig A.l. 
in lRamstedt et al.Ll2008l) . 

Since (1) our data set holds observations up to the high ex- 
citation 7 = 6-5 and 7 = 7-6 rotational line transitions 
of CO, and we want (2) to include the shape of the line pro- 
files, represented by parameter /3, in the estimates of M, and 
(3) to extend the computations to higher mass-loss rate values 
(M> 1 X lO^^Mgyr-'), we elaborated on the studies discussed 
above by performing a grid computation covering more than 
15 000 models. 



3.3. Description of ttie grid 

To study the influence of different stellar and envelope param- 
eters — like stellar temperature, stellar radius (or stellar lumi- 
nosity) and dust condensation radius — on the CO line intensi- 
ties, the input parameters were varied to cover the relevant pa- 
rameter space for evolved low-mass targets (see Table |2]l. The 
values for the /h.co w ere chosen to be 1 x lO''* for O-rich stars 
(Kahane & JuraHl994l) . 3x10 '* for S-type stars (Ramstedt et aU 
[2OO6), and 5x10"'* for C-rich stars dZuckerm an & Dyck. 19861). 
All computations were performed for a 15 m class telescope, 
with beam sizes of 10", 20" and 30" as grid input values. This 
grid spans 15 309 star models, for which each rotational line in- 
tensity was computed for three values of the beamwidth. 

As was done for the observational line profiles, the theoret- 
ical line profiles were fitted with a soft parabola. In the case of 
the theoretical line profiles, however, an extra offset Tf^ should 
be added to the equation, resulting in 



= n 



1 - 



"LSR 







for |v' 
for - vi 



(7) 



10^ 



102 




d = 300 pc 
d = 1000 pc 
d=3000 pc 



10-2 10- 



10° 



IQi '02 10^ 
[K km/s] 



1 0= 



Fig. 5. The continuum corrected integrated intensities 7-^ '^'^ pre- 
dicted with GASTRoNOoM are compared to the values 7*^ cal- 
culated with Eq. [8] Values for the exponents 5, y in Eq. [8] are 
specified in the upper left corner. The full line represents equal- 
ity of 7^''^' and 7°, the dotted lines show a factor 3 difference 
w.r.t. this relation. For all three distance values assumed in the 
model grid, the estimates are plotted in a different color: black 
for d = 300 pc, red for d = lOOOpc, and blue for d = 3000 pc. 
See Sect.[34]for further details. 



where T is the line intensityQ in main-beam temperature. The su- 
perscript 'g' designates that we are dealing with theoretical pro- 
files that were computed using the GASTRoNOoM-code. The 
reason for the offset is that, in case of optically thin envelopes, 
one still detects the almost unattenuated stellar continuum (see 
Fig. |4]i. This offset parameter is not present in Eq. [T]that was 
used to fit the data. During the data-reduction process, the offset 
is removed from the spectra when baselines are subtracted, since 
the observer does not know if and how much stellar light pene- 
trates the envelope. When confronting observational data with 
theoretical integrated line intensities, one therefore has to com- 
pare 7* - 2vc„T^ to 7, where 7 is the integrated line intensity 
jT(v)dv in Kkms"'. The quantity 7' - 2vooT^^ will from now 
on be referred to as 7' '^^, with the superscript 'cc' denoting the 
continuum correction. Optically thin lines allow one to trace the 
wind acceleration zone of which the emission is visible as a su- 
perimposed small peak around the line centre (see Fig.|4]i. 

The soft-parabola fit is not ideal to reproduce optically thin 
resolved emission lines (see upper panels of Fig. [Hi. The /?- value 
was determined in a way that the soft parabola should yield a 
good representation of the line depression. The integrated inten- 
sity of the fit was not required to match the integrated intensity 
of the GASTRoNOoM-predictions. This method provides infor- 
mation about the degree to which the line profiles are resolved 
by the telescope beam. 

The values in Table |3] show that the dependence of the in- 
tegrated intensity on the distance is strongest. For the unsatu- 
rated regime, the integrated intensity is quite sensitive on the 
terminal velocity and mass-loss rate, while this dependence is, 
as expected, much lower in the saturated regime. The intensity 
is moderately dependent on /3, the stellar temperature, stellar ra- 
dius, and dust condensation radius. 



LSR 



> v° 



'* Throughout the paper, we denote line intensity, expressed as a tem- 
perature, with T, and (velocity-)integrated line intensity with /. 
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"CO 1-0 



(5 = 0.06 
/i=0.76 
a = 0.1 6 
tt{3)= '.69% 
7t(10)= 0.09% 



(5 = 0.06 
/i=1 .27 
a = 0.26 
7t(3)= 5.16% 
7t(10)= 0.44% 



10" 10"- 



"CO 3-2 



(5 = 0.02 
/2=0.96 
a = 0.1 2 
tt{3)= 0.46% 
7t(10)= 0.04 



10"^ r 



(5 = 0.06 
/2=1 .08 
a = 0.29 
7t(3)= 9.67% 
7t(10)= 0.09% 



10"° 10"- 
Mdot^ [^©/yt 



-1 1 

log{Mdot„/Mdot„p„) 



-1 1 

log{Mdot„/Mdot„p„) 



'"CO 2-1 


1 1 

(5 = 0.01 
- /i=0.58 <M 
a = 0.09 

tt(3)= 0.34% ^ 
7t(10)= 0.06% ^ 


' /' ' ' 






- /^■eo 7(3) 

/ / 7t(10) 

. .■■ / 1 


> 

(5 = 0.07 

/i=0.86 

a = 0.28 
= 7.67% - 
= 0.05% ; 

1 < 


"" 10"' 10""^ 10"^ 
Mdot„ [Ms/yr] 


10"' 10 


'"CO 4-3 




1 1 1 ji ' " 

(5 = 0.01 .'' / 

- /i=o.57 <mtl^^Kr - 

a = 0.09 / / 
tt(3)= 0.04% y,mi^^^ 


7t(10)= 0.04% 




'■' . .■■ / 1 1 


(5 = 0.05 
/i=1.04 
a = 0.27 
= 6.16% - 

= 0.08% ; 



10"° 10"- 
Mdot„ [Ms/yr; 



-1 1 

log{Mdot„/Mdot„p„,) 



log{Mdot„/Mdot„p„,) 



"CO 6-5 



= 0.02 
= 0.66 
= 0.12 

3)= 0.28% 
10)= 0.05% 




(5 = 0.04 
u,=0.94 
a = 0.23 
7t(3)= 3.08% 
7t(10)= 0.04% 



lO"" 10"= 
Mdot„ [Mo/yr] 



log(Mdot„/Mdot,„„„,) 



XO 7-6 



(5 = 0.02 
u=0.55 
a = 0.1 3 
tt(3)= 0.64% 
tt(10)= 0.04% 



7t(3 



(5 = 0.04 
u=0.65 
CT = 0.22 
1 .83% 



7t(10)= 0.04% 



10"° 10"= 
Mdot„ [Mo/yr] 



-1 1 

log{Mdot„/Mdot,„„„,) 



Fig. 6. Comparison between the input value Minput and M-q resulting from using Eq.|9]for different rotational CO line transitions. 
Estimates for the non-saturated regime are plotted in black, for the saturated regime in grey. The dotted line represents equality 
of the input and estimated mass-loss rates, the dashed lines show a factor 3 difference w.rt. this relation. The listed 6, cr, and /i 
are the mean, the standard deviations, and the maximum absolute values of log(MD/Mnput), respectively. 7r(3) and 7r(10) show the 
percentage of grid models that produce values of Mq deviating more than a factor 3, resp. 10, from Minput. The histogram on the 
right of each panel shows the peak-normalised distribution of log(MD/Minput)- The vertical dashed lines in these histograms again 
show a factor three difference between Mq and Minput- The black and grey histograms represent estimates for the non-saturated, and 
saturated regimes, respectively. 
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Table 3. Values for the exponents i, y in Eq.|9]for each '^CO (J ^ J -1) transition. The boundary values for saturation (see text for 
further details) are given in the third column. The i, y are calibrated such that the mass-loss rate Mu is in units of yr"'. Standard 
deviations cr, maximum absolute values fi, and mean values 6 of log(MD/Minput) are given in the last three columns. 



Regime 


Transition 


Imb X 
[Kkms"' xpc^] 


s 


i.j 






S4.J 


^5J 








Sgj 


'5|0.J 


cr 




6 


unsat. 


1 _ 





< 1 X 10' 


9.73 X 10^ 


16 


1.25 


0.31 


-1.78 


-1.29 


1.03 


0.33 


0.10 


0.21 


-0.01 


0.16 


0.76 


0.06 


sat. 


1 _ 





> 1 X 10' 


9.16 X 10 


-9 


0.78 


0.50 


-1.22 


-0.34 


0.44 


-0.66 


0.28 


0.37 


-0.03 


0.26 


1.27 


0.06 


unsat. 


2 - 


1 


< 2 X 10' 


1.00 X 10- 


15 


1.26 


0.05 


-1.96 


-1.26 


0.79 


0.05 


0.24 


0.10 


0.07 


0.09 


0.58 


0.01 


sat. 


2- 


1 


> 2 X 10' 


9.04 X 10^ 


10 


0.75 


0.43 


-1.47 


-0.08 


0.39 


-1.00 


0.35 


0.41 


-0.00 


0.28 


0.86 


0.07 


unsat. 


3- 


2 


< 9 X 10' 


9.93 X 10^ 


14 


1.14 


0.13 


-1.92 


-1.39 


0.78 


-0.25 


0.40 


0.29 


-0.11 


0.12 


0.96 


0.02 


sat. 


3- 


2 


> 9 X 10' 


9.34 X 10- 


10 


0.61 


0.31 


-1.67 


-0.00 


0.36 


-1.00 


0.37 


0.43 


0.01 


0.29 


1.08 


0.06 


unsat. 


4- 


3 


< 9 X 10' 


9.69 X 10- 


15 


1.17 


0.04 


-1.98 


-1.26 


0.79 


-0.38 


0.39 


0.32 


0.02 


0.09 


0.57 


0.01 


sat. 


4- 


3 


> 9 X 10' 


9.30 X 10- 


10 


0.63 


0.25 


-1.77 


-0.06 


0.38 


-1.00 


0.42 


0.45 


0.03 


0.27 


1.04 


0.05 


unsat. 


6- 


5 


< 1 X 10** 


9.43 X 10- 


14 


1.12 


0.01 


-1.99 


-1.19 


0.75 


-0.48 


0.46 


0.42 


-0.04 


0.12 


0.66 


0.02 


sat. 


6- 


5 


> 1 X 10** 


9.52 X 10- 


10 


0.69 


0.15 


-1.86 


-0.22 


0.49 


-1.00 


0.52 


0.46 


-0.11 


0.23 


0.94 


0.04 


unsat. 


7 - 


6 


< 2 X 10** 


9.44 X 10- 


12 


1.09 


0.01 


-1.99 


-1.06 


0.77 


-0.42 


0.57 


0.52 


-0.02 


0.13 


0.55 


0.02 


sat. 


7 - 


6 


> 2 X 10** 


9.44 X 10- 


10 


0.64 


0.12 


-1.87 


-0.21 


0.47 


-1.00 


0.51 


0.47 


-0.15 


0.22 


0.65 


0.04 


unsat. 


10- 


-9 


< 2 X 10** 


9.73 X 10- 


10 


0.98 


0.00 


-2.00 


-0.72 


0.68 


-0.47 


0.63 


0.60 


0.01 


0.13 


0.45 


0.02 


sat. 


10- 


-9 


> 2 X 10** 


9.41 X 10- 


10 


0.58 


0.03 


-1.93 


-0.08 


0.49 


-1.00 


0.51 


0.47 


-0.15 


0.24 


0.67 


0.05 


unsat. 


16- 


15 


< 2 X 10** 


9.84 X 10- 


10 


0.87 


0.01 


-1.98 


-0.57 


0.61 


-0.85 


0.65 


0.61 


0.02 


0.16 


0.52 


0.03 


sat. 


16- 


15 


> 2 X 10** 


9.20 X 10- 


10 


0.48 


0.02 


-1.92 


0.00 


0.47 


-1.00 


0.46 


0.43 


-0.09 


0.31 


0.92 


0.10 



3.4. Grid results 

To estimate the minimisation function to the integrated intensi- 
ties, the dependence of the integrated intensity F on the different 
input parameters was inspected, and a power-law relation be- 
tween the integrated intensity and all parameters was assumed. 
To estimate the outcome variable we therefore have min- 
imised the function 



CO 



(8) 



The minimisation was carried out for each '^CO Une transition , 
using the generalized reduced gradient method d Lasdon et al.L 
119781) . The coiTelation between the GASTRoNOoM output vari- 
able F''^'^ and the integrated intensity derived from Eq. [8] is 
not bad, but a clear saturation effect is visible for all transitions 
when F''^'^ exceeds a certain boundary value (see Fig. |5]l. These 
boundary values are higher for higher J. A higher mass-loss rate 
not only leads to an increase of the cooling of the circumstellar 
gas, but also to a d ecrease in drift velocity between the dust and 
the gas ( Decin et al.. 2006). This leads to a weaker dependence 
of the integrated intensity of the CO lines on M, visible as sat- 
uration. Not accounting for this saturation effect typically leads 
to underestimating M by a factor of 3 - 10 for M > 3 x 10"^ Mq 
yr-'. 

Moreover, it is visible from Fig. |5] that the boundary values 
for saturation are distance dependent, with lower boundary val- 
ues for larger distances. This dependence was removed by min- 
imising Eq. [8] again, but with both the left-hand and the right- 
hand side multiplied with a factor d^. The coefficients Sij for 
7 = 1 - 0, 2 - 1, 3 - 2, 4 - 3, 6 - 5, 7 - 6, 10-9, 16 - 15 
and the boundary values of Imb x for saturation are listed in 
Table [31 

Replacing 7*^ in Eq.[8]with the observed integrated line inten- 
sity /mb yields an estimate for the mass-loss rate, Mu, provided 
that the other parameters are known: 



(9) 



In Fig. |6]we show a comparison between the mass-loss rate 
A^input used as input for the grid models, and the mass-loss rate 
Md, calculated from the CO line intensities that result from the 
models, using Eq. |9] In Table [3] and Fig. |6]we listed the mean 
values 6, standard deviations cr, and maximum absolute values 
jj. of the quantity \og{Mul Mia^^d- All listed <5-values are posi- 
tive, meaning that on average we slightly overestimate the input- 
mass-loss rate. This is however limited to overestimates of a few 
percents only, except for the cases of 1 - (15 %, 5=0.06), and 
2 - 1 in the saturated regime (17 %, 5=0.07). These deviations 
are still well within the standard deviation cr and therefore not 
significant. In the panels of Fig.|6]we also mention the quantities 
7r(3) and ;r(10). These are the percentages of models in the grid 
that result in a value of Md that deviates by more than a factor 
3, resp. 10, from the input mass-loss rate. The cr-values indicate 
uncertainties of a factor ~1.3 for the unsaturated regime, and 
a factor ~1.9 for the saturated regime when comparing to 
Miapui- When applying Eq.|9]to data, the error bars on the mea- 
sured parameters, e.g. T*, add to these uncertainties, increasing 
the typical uncertainty to a factor of three. 



The spread on the estimates reflects the interplay between the 
different stellar parameters and the line intensities, and the large 
parameter space covered by the grid. In the saturated regime 
the spread is somewhat larger than in the unsaturated regime, 
with standard deviations cr ranging up to 0.16 and 0.31 in the 
unsaturated and saturated regimes, respectively. These values 
are still rather low, and are also reflected in the values of n{3) 
and 7r(10) that do not exceed 10% and 0.45% in the satu- 
rated regime, clearly showing the very small number of outliers 
present in our results. We point out that we found no obvious 
relations between the (few) grid points showing the larger val- 
ues of I log(MD/Minput)|. Therefore, we do not attach larger in- 
trinsic uncertainties to specific input parameters. The quantities 
cr, ju, 6, 71 indicate a high accuracy of the estimator for the com- 
plete range of mass-loss rates included in the grid (Table|2]i. 
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Table 4. Standard deviations cr and maximum deviations /i of 
log(MD/Minput) for the different numbers of parameters included 
in the formalism. Values for both the unsaturated and saturated 
regimes for '^CO (7 = 2 - 1) are shown. 



'-CO (7 = 2-1) 


cr 










unsat. 


sat. 


unsat. 


sat. 


all parameters 


0.09 


0.28 


0.58 


0.86 


no /3, i?,„„e, 


0.10 


0.42 


0.62 


1.13 


nOyS, i?i„ner, T^,,R^ 


0.19 


0.46 


0.78 


1.32 



3.5. Discussion on tine derived mass-ioss rate formuiae 

3.5.1 . Varying the number of free parameters 

As was done in other studies (e.g. lRamstedt et ani2008h . one can 
lower the number of free parameters in Eq. |8] This is illustrated 
in Fig.|2]for the 7 = 2-1 transition. We show the estimates in 
two cases, omitting the dependence of the integrated intensity on 
y6 and Riana- (left panel) and on /3, Rinmi, Ti,, and (right panel). 
Mainly for the mass-loss rates exceeding 10"^ Mq yr"', the esti- 
mates exhibit larger uncertainties. An overview of the values for 
the standard deviations cr, and maximum absolute values jj of 
log(MD/A^f input) is given in Table |4] for the different cases. Both 
fi and cr increase significantly when lowering the number of pa- 
rameters. This implies that it is better to include a well-founded 
estimate for some parameter, like /3 or /dinner, than to omit it from 
Eq.El 

3.5.2. Sensitivity analysis 

To further assess the quality of our estimator, we focus on the 
influence of(l) input parameters for Eq.|9] such as effective tem- 
perature, luminosity, dust condensation radius, photospheric CO 
abundance, and distance, (2) the envelope's outer radius, and (3) 
the gas kinetic temperature. 

Input parameters: We tested the influence of d, T*, (R*), 
/h,co and dinner ou the M-cstimatcs for 7=1-0,2-1,3- 
2, 4 - 3, 6 - 5, 7 - 6. We assumed a standard stellar model with 
Minput = 10-^ Mo yr-i, = 2000 K, R^ = 3.14 x 10'^^ cm, 
■dinner = 3/?*, /h,co = 10""*, d - 300pc, and a fixed beamwidth 
Ob = 10". The integrated intensities /mb are those calculated 
with the GASTRoNOoM-code for this model. In Fig.[8]we show 
the influence of the different parameters on the M-estimates by 
varying d, /h,co, T*, L* and Rinmr one by one while keeping all 
others constant at the model values. We point out that a slight 
gradient is present in the M-estimates for the standard model, 
i.e. for higher 7 we seem to estimate somewhat lower values of 
M. The deviations and the spread of the estimates are, however, 
well within a factor three from the input mass-loss rate. Since the 
latter is used as a criterion to decide upon constancy of M, we 
can state that the estimator can very well reproduce the input-M. 

As expected, the distance has the strongest influence on 
the estimated M-values, with a similar increase or decrease in 
mass-loss rate for all transitions. Altering fn.co has a pure scal- 
ing effect on the estimates, with a lower M for a higher /h,co- 
Changing T*, L*, or /?innei leads to the largest differences for the 
higher-7 transitions, since changing these parameters affects the 
layers closest to the star more strongly than it does the more out- 
ward layers. This effect is reflected in the larger absolute values 
of exponents sjj, ig y and S()j in Eq.|9] and by the fact that the 
mentioned gradient is slightly shallower or steeper when con- 



sidering other values for T*, L*, or /^innei- This implies that the 
estimates based on the higher-7 transitions suffer from larger un- 
certainties and that more value can be attached to the estimates 
based on the lower-7 transitions. 

Inspecting the different panels in Fig. [8] we find no indica- 
tions for changes in the M-estimates with an opposite character 
for the different transitions, in the sense that, e.g., the estimates 
for some 7 would go lower, while those for other transitions go 
higher. For a constant mass-loss rate we will therefore only see a 
gradient (if any) and not a random distribution of the estimates. 
This implies that a faulty parameter assumption will only lead to 
a general shift of the estimates and is unlikely to change the in- 
terpretation of possible variability or constancy of the mass-loss 
rate. 

Outer envelope radius: The influence of the outer radius of the 
CSE, Router, can not be tested in the same way since it is not an 
input parameter in Eq. |9] This variable has therefore no direct 
influence on the estimates. It did however influence the deriva- 
tion of the coefficients s/j, since the geometrical extent of the 
envelope affects the integrated intensities and the shapes of the 
emission lines. The transition most sensitive to changes in Router 
is 7 = 1-0, implying that it is a safer choice to use M-estimates 
based on the 7 = 2-1 transition. 

Clumping of CSE material causes a larger 7?outer than a 
smooth mass outflow because of shielding properties of the 
clumps. This implies that higher integrated intensities can be 
reached with the same mass-loss rate. Therefore, if clumping is 
present, we could be overestimating the mass-loss rate. This will 
again be most strongly reflected in the 7=1-0 estimates. 

Gas kinetic temperature: In Fig.|9]we compare gas kinetic tem- 
peratures T, which have been consistently computed with the 
GASTRoNOoM-code, and integrated intensities F "^^ of transi- 
tions 7 = 1-0 up to 7 = 7-6 for two different values of 
Mnput (IO^^Mq yr"' and 10"^ Mq yr"') and five different ef- 
fects in cooling. All other model parameters are the same as 
those mentioned above for the standard model. The default cool- 
ing scheme, which was used to derive the formalism presented 
in this paper, includes among others cooling due to rotational ex- 
citation of H2O. The cooling schemes we considered (1) ignore 
the H20-cooling, (2) include H20-cooling (default), (3) include 
H20-cooling ten times stronger than assumed by default, (4) as- 
sume a power law for the gas kinetic temperature 

T(r)^TARJry (10) 

with e = 0.5, and (5) assume a power law with e = 0.7. 

We point out that, except for the cases in which a power law 
was used to describe T(r), the temperature can increase or de- 
crease with increasing radius r, i.e., instead of cooling there can 
be heating of the envelope; see Fig.|9l 

We focus on the H2O rotational cooling since this part of the 
cooling is still subject to large uncertainties. In GASTRoNOoM, 
H2O is modelled as a three-level system, which is a gross sim- 
plification of the actual structure of the molecule. The rota- 
tional rate c oeffic ients used here are based on the re s ults o f 
Green et al.' ('1993'), as are those of lNeufeld & Kaufmanl(ll993h . 
Faure & Josselin (2008) recently c alculated new rotational rate s 
and compared them to those of iNeufeld & KaufmanI (Il993h . 
finding higher cooling rates for temperatures < 1500 K and 
lower cooling rates for temperatures > 1500 K. They also found 
that the rotational rate coefficients can be uncertain by factors 
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Fig. 7. Comparison between the input value Minput (ordinate) and Md resulting from using a formula analogous to Eq.|9] (abscissa) 
for different rotational CO line transitions, but omitting the dependence of the mass-loss rate on Rinnei and p (left panel), and on 
dinner, j8, , and (right panel). Estimates for the non-saturated regime are plotted in black, for the saturated regime in grey. 
6, a, fi, 7r(3) and ;7r(10) are as in Fig|6l The dotted line represents equality of the input and estimated mass-loss rates, the dashed 
lines show a factor 3 difference w.r.t. this relation. The histogram on the right of each panel shows the peak-normalised distribution 
of log(MD/A^input)- The vertical dashed lines in these histograms again show a factor three difference between Md and Minput- 
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Fig. 8. Sensitivity check of the used formalism for five parameters: d, /h,co, T*, (R*), and Rinmi-. The upper left panel gives the 
estimates based on the input parameters of the model and the output values of the integrated intensities 7^ '^'^. The grid values for the 
different parameters and the colour coding are specified in the legend of every panel. The dash-dotted line in all panels marks the 
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(a) Gas kinetic temperature T versus radius r for M = 10 ' Mq yr ' (full 
line), and 10"^ Mq yr"' (dotted line). 
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(b) Integrated intensities F-'^'^ for transitions 7 = 1-0 up to 7 = 7- 6. 

Fig. 9. Gas kinetic temperatures (a) and integrated intensities 
(b) for different Minput-values and different cooling mechanisms. 
The colour coding in both panels is as follows: black was used 
for omission of H20-cooling, blue for inclusion of H20-cooling, 
green for inclusion of H20-cooling ten times stronger than by 
default, yellow for a temperature power law, as in Eq. [TOl with 
e = 0.5, and red for e - 0.7. In the lower panel, the different val- 
ues of the input-mass-loss rate are indicated with diamonds and 



triangles for M =10 Mq yr 
See discussion in Sect. 13.5.21 



and 10 Mq yr , respectively. 



of a few up to an order of magnitude, which is reflected in the 
uncertainties on the cooling rates. 

The contribution of H20-cooling is most important in the in- 
ner wind regions, since farther out in the envelope CO-rotational 
cooling and especially cooling due to adiabatic expansion dom- 
inate. This implies that the gas kinetic temperature in the inner 
wind regions is more uncertain than in the rest of the envelope. 

As is visible in Fig. |9(b)| decreases for increasing 
H20-cooling. Again, this effect is more explicit for higher J. 
Overestimating the H20-cooling therefore leads to an underesti- 
mate of the mass-loss rate, considering the relation between the 
integrated intensities and the M-estimates. 

Since the excitation regions of the higher-7 transitions are 
partly situated in the inner wind regions, the integrated inten- 



sities of these transitions will be subject to larger uncertainties 
than those of the low-7 transitions 7=1-0 and 2-1. This 
different degree of influence on the integrated intensities for the 
different 7-values is very clear in Fig. |9(b)[ where the spread on 
the integrated intensities of the transitions 7 — > 7 - 1 increases 
for higher 7. A second implication is that the mentioned gradient 
of (slightly) lower M for higher 7-values is possibly linked to the 
uncertainties on the H20-cooling in the inner regions. If this is 
indeed the case, then the H20-cooling rates used in our models 
are overestimates of the actual rates. 

Considering the current models, the uncertainties on the in- 
tensities of higher-7 transitions (as discussed above), and taking 
into account that the intensities of the 7=1-0 transition can be 
affected by e.g. masering and clumping, we put forward that the 
7 = 2-1 transition is likely the most reliable transition to use in 
estimating mass-loss rates. 

3.5.3. Comparison to lRamstedt et al.l <2008|) 

In benchmarking our formalism against mass-loss rate estima- 
tors presented in the literature, we focus on the one derived by 
Ramstedt et al. (2008). A di rect compar i son to the estimators of 
.Knapp & Morris, (.19851) or iLoup et all (Il993h is not possible, 
since these have the main-beam temperature at the line centre, 
Tmb.c, as a fundamental parameter, while we consider the inte- 
grated intensity /mb and the shape of the lines. 

In Fig.[TO]we show the mass-loss rate estimates derived for 
transitions 7 = 1-0, 2-1, 3-2 and 4-3 via Eq. |6] us- 
ing the same parameter grid as was used to de rive the formalism 
presented in this paper. Ram stedt et al.l (1200 8) mention that their 
estimator Mr is only valid for mass-loss rates between 10"^ and 
10"^ Mq yr so we ma rked those outside this interval in grey, 
for the sake of the claritv. lRamstedt et al.l (l2008l) also note the re- 
quirement of unresolved envelopes to use their estimator to high 
accuracy, and we therefore eliminated the grid points that result 
in resolved envelopes by determining the CO-photodissociation 
radius, Rco, and su bsequently the oute r rad ius, isomer, of the en- 
velope according to lStanek et al.l(ll995h and lSchoier & OlofssonI 
(,200L) : 



Rco = 5.4 X 10' 



M 

To 



-6 



0.65 



,15 ./h,co(-R*) 

X X 7- 

v«, 4x10-4 



0.55 



.7.5X10-(:J|) 



„ = 2.79x1 Axi^ 



cm 

0.09 



^lo-f- 

/h,co('-) = /h,co(^*) X exp 



(11) 

(12) 
(13) 



The CO-photodissociation radius and the outer radius of the en- 
velope are the radii at which the CO abundance has dropped to, 
respectively, 50 % and 1 % of its photospheric value, /h,co(^*)- 
a describes the abundance profile throughout the envelope. We 
adopted the Minput-value to calculate Rco and /?outer- When the 
angular size of the envelope with radius /?outei exceeded the beam 
size 6b, we excluded the respective grid point. 

The mean values 6, standard deviations cr, and maximum ab- 
solute values /i of log(MR/Minput) are printed in each panel of 
Fig. [To] both for the full grid and for only those grid points in- 
side the interval [10^^; IO^^JMq yr"'. Very high values of yu 
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Fig. 10. Minput versus Mr calculated from the grid also used to construct Fig. |6] and Eq. |6] The dotted line indicates where the 
quantities are equal, the dashed lines show a factor three difference w.r.t. this relation. The listed 6, fi, and cr are, respectively, the 
mean values, the standard deviations, and the maximum absolute values of log(MR/Minput)- We specify these quantities for the full 
grid (plotted in grey), and for the grid points inside the interval [10^^; 10"^] yr ' (plotted in black). The histogram in each panel 
shows the peak-normalised distribution of log(MR /Minput)- 



are reached when the complete grid is considered. Confining the 
statistics to the interval [10^^; 10"^] Mq yr"', we find that these 
values are much lower. For all transitions fi and cr are higher 
than the values mentioned in Table [3] derived for our estimator. 
There is also a significant shift to higher M- values, for all but the 
/ = 1 - transition, using the estimator of' Ramstedt et al.l (l2008l) 
in Eq.|6] This shift is indicated by 5-values significantly deviat- 
ing from zero (6 > cr) and is expected due to thedifferencesjn 
the tem perature structure adopted by us and by iRamstedt et al.l 
(l2008h . The uncertainties on the cooling, especially by rotational 
excitation of H2O, are still rather large, especially so in the inner 
wind region, as discussed in Sect. 13.5.21 Furthermore, our ap- 
proach utilises a higher number of parameters and expands the 
validity range of the M-estimator by correcting for saturation of 
the CO lines, ensuring that we get small intrinsic uncertainties 
for our estimator. 



4. Stellar parameters 

4. 1 . Temperature, luminosity and distance 

Since we want to obtain empirical relations between wind pa- 
rameters (e.g. M, v'co), molecular-line parameters (e.g. line inten- 
sity), and basic stellar parameters (e.g. T*, L*), the latter need 
to be well determined. To draw meaningful conclusions from 
sample statistics, a highly homogeneous determination of basic 
stellar parameters is needed. Effective temperatures were derived 
from the dereddened V-K colour, and luminosities from period- 
luminosity relations. A description of the methods used to ob- 
tain these parameters is given in Appendix|C] The results for the 
sample are given in Table|6] together with the evolutionary stage, 
chemical type and pulsational type of the objects. 

4.2. Inner envelope radius 

As discussed in Sect. [3] the inner radius of the circumstellar en- 
velope, ./dinner, IS uscd in determining the mass-loss rate. In es- 
timating M for the sample targets. Thinner is calculated assuming 
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Table 5. Overview of sample targets for which no mass-loss rate 
estimates were made. The third column lists the evolutionary 
stage of each target. 



IRAS 


Other identifier 


Evolutionary stage 


r\/^ 1 in^/^ 
Uol /o-lUJo 


Red Rectangle 


r-ALrr> 


\j / jyy- I'^jj 


Calabash Nebula 




lUiy /-J /jy) 


Oc>C UooUo-UUjUy 


r-ALrD 










Alpha Sco 




1 / 1 jU-OzZ4 


Cotton Candy Nebula 


D A r^R 


MAA'X 0Q/1Q 




D Ar^R 
r-ALjD 


1 / jui-zdjo 


V'4-ZUi ogr 


AvjJj 




Gomez Nebula 


I oU 


18100-1915 


OH 11.52 -0.58 


AGB 


18257-1000 


V441 Set 


AGB 


18308-0503 


AFGL 5502 


YSO 


18327-0715 


OH 24.69 +0.24 


AGB 


18361-0647 


OH 25.50-0.29 


AGB 


18432-0149 


VI 360 Aql 


AGB 


18460-0254 


VI 362 Aql 


AGB 


18488-0107 


VI 363 Aql 


AGB 


18498-0017 


VI 365 Aql 


AGB 


19067+0811 


VI 368 Aql 


AGB 


191 10-H 1045 


KJK G45.07 


Hll-region 



that the dust temperature at this radius is equal to the conden- 
sation temperature, and that 



Yr) ^^^^ 



,2/(4+s) 



with s K 1 (Olofsson in iHabing & Olofssonl l2003h . The con- 
densation temperature is taken to be 1500K for O-rich and S- 
type stars and 1 200 K for C-rich stars. 



Table 6. Ov erview of the evolu tiona ry stage, atmospheric chem istr y, pulsational type, spec tral type, and ste l lar pa rameters of the sample stars. Spectra l types were taken from (i) the SIMBAD 
database, (iP lKwok et"al] | |1997|) . (iii) lFeast & Whitelockl ilOM) . (iv) lJusttanont et all i ll 9961) . (v) ISuarez elal] ( 120061) . and (vi) ICastro-Carrizo et all ( 120071) . For each star the following fundamental 
parameters are specified: distance d, period of pulsation P, effective temperature T^g, luminosity Lt,, and stellar radius R^,. T^ff and were determined according to the methods presented in 
Appendix|C] 

References to the values adop ted from the literature are: (l)lMassey et al1(l2006l). (2)ISch6ier et al](l2007h. (3)lLevesaue et al.l (l2007l). (4) |Matsuura et al.| ( |2002|) , (5) lpecin et"al] (l2008Natureexlabbl) . 
(6)i.Tevssier et al. (2006), (7) IWhitelock & Catchpolj ( Il985h . (8) lLevesque et al.l ( l2005h . (9) ISch6ier et alj (l2005h . ('lO)lFaiindez et al. (2004), (11) de Jaeei (1998*). (12) lLikkel et al.l (Il99lh . (13) 
iGroenewegen et al.. (,19991) . (V - Kh is the V — K colour that was dereddened according to the method discussed in Appendix |B] A',^ean is the mean K-band magnitude. 
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Table 6. Continued 
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Table 7. Mean mass-loss rates M derived with Eq. |9]for all AGB, Post-AGB, supergiant, and liypergiant stars from the sample. Targets with line profiles that show clear deviations from the 
soft-parabola profile are marked with (!). a(M) indicates the spread of the M-estimates if a minimum of three diagnostic lines is available. If variability is detected in the estimates, this is indicated 
in the fourth column. If no variability could be seen in the estimates, this is indicated with Constant. The fifth column lists some general remarks. Mines is the number of lines used. The telescopes 
with which the line data were obtained are listed in the seventh column. Romki is given in the last column. 
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Fig. 11. M-estimates from Eq.|9]for the AGB stars in the sample. The dashed line indicates the mean value M listed in Table |7] the 
dashed line indicates the mean of the estimates. The shaded area indicates the typical uncertainty region around the mean value, 
which is a factor of three. 
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Fig.[TTl(continued) 
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5. Results and discussion 

5.1. Mass-loss rate M 

For all targets in the sample for which we could fit at least one 
rotational line profile with a soft parabola (Eq. [1]) we listed esti- 
mates for M in Table|71 The estimates from all '^CO lines for the 
sample AGB stars are shown in Fig.[TT] For the supergiants and 
hypergiants they are shown in Fig. [12] In case the soft-parabola 
fit procedure yielded a negative fi, the exponent siq j in Eq.|9]was 
assumed to be zero. 

The listed estimates are averages over all values obtained 

via Eq. |9] for the fitted line transitions, and are written as M. 
In Table|7] we also listed the number of lines used and the spread 
of the estimates around the mean. 

We did not attempt to fit the line profiles for U Ant and S Set, 
since these are composed of emission from the inner part of 
the CSE and from the detached shell at larger outflow veloci- 
ties, respectively. As discussed in Sect. 12.3.21 some other targets 
show clear deviations from the fitted soft-parabola profile. The 
estimates for these targets, e.g. EP Aqr, are therefore subject to 
larger uncertainties than the spreads listed in Table |7] and have to 
be interpreted with caution. 

If the mass-loss-rate estimates of all available lines are 
within or close to the typical factor three uncertainty region, it 
is reasonable to state that the mass-loss rate of the target has 
been constant throughout the regions of the envelope sampled 
by these transitions. For most AGB stars in the sample this is 
indeed the case, but there are some exceptions indicating a pos- 
sible variability in mass loss in the part of the envelope that is 
traced. 

5.1.1. Variability 

Criterion: Since higher-7 transitions trace deeper, more recently 
formed envelope layers, a criterion to decide upon variability 
should consider the behaviour of the M-estimates in view of the 
7-values. The presence of a picket-fence structure in the esti- 
mates, i.e. with M varying strongly between transitions could 
suggest variability of the mass-loss rate, but on rather short time 
scales of some hundreds of years. The presence of a strong over- 
all gradient can also imply mass-loss-rate variability: e.g. de- 
creasing M for increasing J can be indicative for a mass-loss 
rate that has decreased over time. When discussing the sensitiv- 
ity of our results to the modelling of the cooling of the envelope 
in Sect. 13.5.21 we pointed out that our estimates might intrinsi- 
cally hold a shallow gradient. The differences associated with 
these model uncertainties between the highest and lowest esti- 
mates (for a constant Minput) are, however, not expected to ex- 
ceed a factor three for J - 1-0 up to 7 = 7- 6. Therefore, if 
a much steeper gradient is found in the estimates for some tar- 
get, with a total coverage exceeding an order of magnitude in M, 
we are therefore inclined to state that the mass-loss rate of the 
respective target has varied over time. 

Variability in the sample: The estimates for V669Cas, shown 
in Fig.[TT] exhibit the picket-fence type of trend and suggest that 
variability in the mass-loss rate is traced for this star 

The estimates for NV Aur and GLMP 1048 based on data of 
the 7 = 2-1,3-2, 4-3, 6-5 transitions, the estimates for 
RSclandIRC-H20370,for7 = 1-0, 2-1, 3-2, 4-3, 7-6, and 
the estimates for CW Leo, for 7 = 1-0, 2-1, 3-2,4-3,6-5 
exhibit a spread of about an order of magnitude. This spread is 
within the uncertainties we expect around a M-value and there- 



fore points to a constant mass-loss rate throughout the traced 
regions of the envelope. The estimates for RW LMi, LP And and 
xCyg with J - 1-0, 2-1,3-2, 4-3, 6-5, respectively 
cover factors about ten, twenty and thirty in M. The estimates for 
RXBoo, for 7 = 1-0, 2-1,3-2, 6-5, cover a factor about 
forty in M. In all these cases, we find a strong downward gra- 
dient, which most likely indicates that the more inward regions, 
i.e. the excitation regions of the higher-7 lines, were produced 
by a lower mass-loss rate than the outer regions. A strong gradi- 
ent is also found for both WX Psc and LL Peg for J - 1 - up 
to 7 = 6 - 5. For both targets, however, the M-estimate based on 
the 7 = 7-6 transition is higher by about an order of magnitude 
than the one based on the 7 = 6-5 transition, indica ting a recent 
increase of the mass-loss rate, as was proposed by iDecin et al] 
('2007') for WX Psc. Qualitatively, this phenomenon is indepen- 
dent of the modelling of the HaO-cooling, i.e. an increase would 
still be present if another cooling mechanism was considered in 
deriving the estimates. Obtaining a result corresponding to de- 
tailed modelling using an estimator as in Eq.|9]can be used as a 
benchmark for the quality of the estimator 

It is likely that the OH/IR stars and the more extreme Mira 
stars in the sample have already undergone mass-loss modu- 
lations, but to verify this more data are needed, especially of 
7 = 1 - 0, 6 - 5 and 7 - 6. 

Data quality: The IRAM measurementflin our sample are sub- 
ject to somewhat larger uncertainties than our own APEX and 
JCMT d ata since som e questions have a risen about the data re- 
duction. iDecin et aP dToOSNatureexlabah already reported large 
uncert ainties on the IRAM data of RHya and ISkinner et al.l 
(fT999l) mentioned significant deviations between different mea- 
surements of the same lines towards CW Leo. The IRAM data 
ISkinner et al.l gathered from the literature and archives seemed to 
suffer from large uncertainties linked to pointing and calibration. 
They also discussed inconsistencies in the literature concern- 
ing the reported absolute line intensities. Very often confusion 
or unclarity arises about e.g. conversion from antenna tempera- 
tures to main-beam-brightness temperatures. For these reasons, 
we considered leaving the IRAM measurements out of the data 
set. In case of LL Peg this leads to a decrease of the estimated 
mass-loss rate with a factor 5.3, for IRC +20370 with a factor 
3.9. For IRC +10420 and R Hya omitting the IRAM data has no 
effect on the M-estimate. For the eight other objects with IRAM 
observations, the decrease is with a factor between 1.3 and 2.7. 
The M-estimates ignoring IRAM data are always lower. 

Omitting the IRAM data has an effect on the gradients seen 
in the estimates and discussed earlier In many cases, the CO- 
sampling is now restricted to the transitions higher than 7=1-0 
and/or 7 = 2-1, making possible trends less significant. In 
case of CWLeo, RXBoo, IRC +20370, RWLMi and LP And 
the supposed trends in the estimates are no longer present. For 
LL Peg the presence of an OSO 7=1-0 measurement still gives 
the trend seen when the IRAM data are included. 

For the oxygen-rich semi -regular variable RX Boo we find 
two estimates that could hint towards variability, i.e. those based 
on the IRAM 7 = 1 - and JCMT 7 = 6-5 measurements. 
The latter is duplicated by a measurement carried out with CSO, 
and the estimate based on that observation is about an order of 
magnitude higher than the one derived from the JCMT. We are 
however inclined to attach more value to o ur JCMT observation , 
since the CSO calibration as described bv lTevssier et alJ (l2006h 

' All IRAM measurements were previously presented in the litera- 
ture. 
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Fig. 12. M-estimates from Eq. |9]for the supergiants and hypergiants in the sample. The dashed line indicates the mean value M 
listed in Table|2l the dashed line indicates the mean of the estimates. The shaded area indicates the typical uncertainty region around 
the mean value, which is a factor of three. 



holds larger uncertainties, such as rather large pointing errors. 
We conclude from this that there is some, however somewhat 
uncertain, downward trend with lower M-estimates for higher 
J, indicating the possibility of a decrease in mass-loss rate over 
time. 

Comparison to the literature: iTevssier et"al] (120061) note that 
a change in the mass-loss rate is needed to reproduce the data 
for three of their sample targets: RHya, RWLMi and ;i'Cyg. 
[Decin et al. (2008Natureexlaba) constructed a mass-loss model 
for R Hya and showed that a detached shell is needed to explain 
the rotation-vibrational lines. The purely rotational lines could 
however be reproduced with a constant mass-loss rate. Our es- 
timates do, indeed, not point to any significant variations in the 
mass-loss rate, as can be seen from Fig. [TT] The M-estimates 
for RW LMi cover one order of magnitude and show a (shallow) 
gradient with lower M for higher /. This could point to a mass- 
loss rate that is presently lower than the o ne that produce d the 
outer layers of the circumstellar envelope. ITevssier et al.! spec- 
ify a difference of only a factor 1.5 between the two mass-loss 
rates, so we can say that our estimates supp ort this hypothe sis of 
mass-loss-rate variability. In case of x Cyg ITevssier et al. l men- 
tion a minimum change of a factor 4 between the former and the 
present mass-loss rate, with the present one being lowest. This 
result is again reflected in our estimates. 

The detached-shell structure of the carbon-rich semi-regu lai- 
variable R Scl was described bv'Bergeat & Chevallied (l2005h as 
the consequence of a former mass-loss rate of ~5.5 X 10"^ 
yr"' and a present mass-loss rate of ~3.5 x 10"^ M© yr"'. The 
range of mass-loss rates covered by our estimates agrees with 
these values. 

The mass-loss rates estimated from the individual CO emis- 
sion lines for WX Psc — shown in Fig. [TT] — are in very good 



agreement with the detailed modelling of the star bv lDecin et al.l 
(2007). They present a profile Mir) exhibiting a high M-value, 
followed in time by a strong decrease and again an increase of 
the mass-loss rate. We do not reproduce the same large differ- 
ence o f three orders in magnitude as mentioned by iDecin et al.l 
(l2007|). A possible explanation is that the excitation regions of 
the individual observed transitions are not restricted to either the 
high-M or the low-M regions. The estimates then represent more 
than one of these regions and are a kind of 'mean' value of M. 
However, the estimates we present in this paper for WX Psc re- 
flect the detailed mass-loss modelling very well. 

IWinters et al.l (I2OOOI) reported on possible time variability of 
the mass-loss rate of CW Leo. The present-day mass-loss rate 
is suggested to be lower than the mass-loss rate that produced 
the outer envelope layers by a factor of up to five. The decrease 
we see in our estimates, however, is by a factor ~30. Also, the 
M- values presented in the literature are about an order of magni- 
tude higher than the M-value we derived. A possible explanation 
for these discrepancies is the H20-cooling uncertainty, leading 
to underestimates of M. Another explanation is that the 7 = 7-6 
transition for the first time reveals a stronger variation in M than 
could be traced with single-dish observations o f CO transitions 
up to J = 4 - 3 or 6 - 5. The model presented bv lRamstedt et al.l 
( I2OO8I) fitting Unes J - 1-0 up to 7 = 4- 3 shows clear de- 
viations from the data. The discrepancies are systematic, i.e. the 
model underestimates the observed intensities more strongly for 
increasing J. This indicates the need for consideration of M- 
variability for this target. 

iNeri et aP (1199 8*) suggested that LL Peg might have a bipo- 
lar detached shell and Mauron & Huggins (2006) reported on a 
spiral structure in the envelope. The latter is indicative for vari- 
ability in the outflow from the central star. 
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Supergiants and hypergiants: The estimates for all five RSGs 
in our sample suggest constant mass-loss rates; see Fig. [12] In 
case of the yellow hypergia nt IRC +104 20, H umphreys et al. 
(|l997[) : lO^ro-Carrizo et all (l2007l) and Ii5jnh-V.-Trung et all 
(|2009|) . among others, have put forward the possibility of a com- 



plex mass-loss history, with several episodes of high and low 
mass-loss rates. This resulted in multiple arc-like structures and 
shells around the central star. For AFGL2343, the second hy- 
pergiant in our sample, we only used the 7 = 3-2 line of '^CO 
to estimate the mass-loss rate, so no concl usions on variability 
in m ass-loss rate can be drawn from this. ICastro-Carrizo et al] 
(1200 7) mention variations in the mass-loss rate for both hyper- 
giants on time scales of ~1000years. In this study, where we 
only consider CO emission lines measured with single-dish tele- 
scopes, we lack sensitivity for the rapid density fluctuations. This 
is due to the overlap of the line forming regions of the different 
rotational transitions. 



5.1 .2. AGB mass loss in view of pulsational type and 
cliemistry type 

The histograms in Fig. [T3] show that the spread on M for AGB 
stars is quite large, but there is a noticeable difference between 
the semi-regular variables on the one hand and the Miras and 
OH/IR targets on the other — see Fig. |13(a)| The spread on M 
for the 16 Miras is rather extended, with values ranging from 
9 X lO-** Mo yr-i (R Leo) up to 6 x 10"^ Ma yr"! (NX Ser) and 
a mean value of 10"^ M© yr"'. 

The mass-loss rates for the 9 OH/IR targets were all esti- 
mated in the range 4xlO"''<M<6x 10"^ Mq yr"', with 
an average of 2 x 10"^ M© yr"'. The estimates for these targets 
seem low, considering that OH/IR stars have optically thick en- 
velopes and can be he avily obscured. However, it has already 
been pointed out by e.g. lDelfosse et alJ (Il997h that the mass-loss 
rates derived from CO lines can be significantly lower than those 
derived from the infrared flux in the case of OH/IR stars. The 
most likely explanation for the discrepancy is found in the recent 
(~1000yr ago) onset of a superwind phase. The latter would not 
be traced by the CO lines available in this data set, except possi- 
bly for WX Psc. 

The SRb targets have lower mass-loss rates, with an av- 
erage of 5 X 10"^ Mq yr"' and the highest value reaching 
2 X 10"^ Mq yr"' (RScl). These results are in accordance with 
the idea that the SR variables are likely progenitors of Miras 
( Whitelock & Feast, 2000; Yesilyaprak & Asian, 2004). 

When considering the histogram representing the AGB tar- 
gets grouped per chemistry type (oxygen-rich, carbon-rich or S- 
type) in Fig. |13(b)| no clear distinction between these groups 
can be made. Though this could be due to the bias of our sam- 
ple towards oxygen-rich stars, there are no clear indications 
for significant discrep ancies in M between the different types. 
iRamstedt et alJ (l2006l) also reported no noticeable differences in 
M-values for the different chemistries. Putting all AGB targets 
together for which M-estimates could be made — irrespective 
of pulsational or chemical type — in Fig. 1 13(c)] we get a mini- 
mum value of 4 X 10"^ Mq yr"' (SRb V438 0ph), a maximum 
value of 6 X 10"^ Mq yr"' (O-rich Mira NX Ser) and a mean of 
10"^ Mo yr"'. 

5.1 .3. AGB mass loss as a function of pulsation period 

The M-values for all AGB stars in the sample — except for 
the RV Tauri star R Set — are plotted versus pulsation period in 



Fig- [M] with error bars indicating the spread as listed in Table]?] 
In case only one or two estimates could be obtained, we assumed 
a spread of a factor three, which is appropriate considering the 
overall spread i n the sample . The red data points in the graph 
were taken from iHe & Cheii] (]200l and cover different types of 
OH/IR stars with a large spread on th e perio ds of pulsation. Five 
objects in the sample of iHe & Chen! (12001 ) are also present in 
our sample. Their pulsation periods and mass-loss rates are in 
accor dance with o ur estimates, except in the case of V669 Cas , 
where iHe&Chenl use P = 1995 days (ISlootmaker et al.Lll985h . 
while we adopted the more recently determined P = 1525 days 
(iGroenewegen et a"l1,ll999l) . In constructing a linear relation be- 
tween the pulsation period and the logarithm of the mass-loss 
rate, we only took into account the data p oints in our s ample — 
i.e. we did not consider the data points of iHe & ChenI ( 1200 lb — 
with periods shorter than 1000 days. The positions in the graph 
of the objects with periods longer than 850 days are represented 
by a constant M of 3.4 x 10"^ Mq yr"'. The linearly increasing 
part of the fit (f < 850 days) was constructed to represent the 
Miras, the short-period OH/IR stars and the AGB semi-regulars. 
The constant fitting value represents a levelling-off or satura- 
tion of the AGB mass-loss rate for stars wit h long periods, as 
proposed by, e.g.. lVassihadis & WoodI (Il993l) . The lack of data 
points for stars with longer periods complicates the determina- 
tion of the cut-off period and the mean M-value for the long-P 
objects. Both were more or less arbitrarily chosen to yield a good 
fit. The fit in Fig.[T4]is represented by 



log(M) = -7.37 + 3.42 x 10"^^ x P 
= -4.46 



P < 850 days 
P > 850 days 



with M in units of Mq yr"' and P in days. Although He & Chenl 
(12001') note that the quality of their derived mass-loss rates may 
be questionable, we find that their data points of both short and 
long-period targets are well represented by our fit. The targets 
that deviate strongest from the fit are semi-regular variables, e.g. 
WHya and UHya, and the Miras NX Ser (O-rich) and WAql 
(S-type). We note that in constructing the fit we did not take into 
account the effect of differences in stellar mass that are likely 
to cause extra sprea d on the mass-loss rate for fixed pulsation 
periods, as shown bv lVassiliadis & Woodl (Il993h . We also point 
out the ve ry good agreement between our fit and the results pre- 
sented by ISchoier & OlofssonI (1200 ih . who produced a similar 
diagram of M versus pulsation period for a sample of 68 op- 
tically bright carbon stars (their Fig. 12), and found a relation 
M cc p2-3 for periods P between 100 and 800 days. 



5.2. '^C/'^C isotope abundance ratio 

For 29 sample stars we have observations of both '^CO and 
'^CO rotational lines, allowing us to estimate the '^C/'-'C 
isotope abundance ratio. This ratio is a tracer o f the star's 



chemical evolutionary state and nucleosynthesis dS opka et al 
1989; Sch"5ier & Olofssonll2000l;lGarcfa-Hernandezet al...2007 



Mil am et al. , 20091) . '^cTs produced during the interpulse phases 
when, after dredge-up, protons are partially mixed in the region 
between the hydrogen and helium burning shells. Reactions of 
the protons with '^C in this intershell region form '^C. In mas- 
sive AGB stars (3 Mq < < 7 Mq), "hot bottom burning" 
(HBB) can be activated when the convective envelope pene- 
trates the hydrogen-burning shell. During HBB, '^C/'^C is low- 
ered, since the net production of '^C is higher than that of '^C 
in the re-activated CN-cycle. For HBB-stars, a ratio '^C/'^C 
around 10 is expected, but could possibly range down to values 
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Table 8. Line intensity ratios /mb('^CO)//mb('^CO) for the sam- 
ple stars give first order estimates of their isotope abundance 
ratio '-C/'-'C. An absolute flux error of 20% for transitions 
7 = 2-1, 3-2, and 4-3 gives an error of 28 % on the re- 
spective line intensity ratios. For 7 = 6-5, the absolute flux 
error is of the order of 30 %, so the error on the line intensity ra- 
tio is about 42 %. The '^C/'^'C-value listed in the second to last 
column is the mean line intensity ratio, corrected with a factor 
0.87 for the difference in line stren gths (Eq.[T5]l . The last column 
(M09) lists the values derived bv lMilam et all i2QQ% from CO 
data. 
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about 3 or even 1 dOhnaka&TsuiiL 119991) . HBB is expected to 
be active in the most massive and luminous oxygen-rich AGB 
stars. In this respect, OH/IR stars — typified by high lumi- 
nosities, large infrared excesses and high mass-loss rates — are 
the best candidates for having very low ^^C/^^C isotope ratios 
dGarcia-Hernandez et aLllToOTl and references therein). 

Obser vational ^^C/^^C -values for AGB stars r ange from as 
low a s 1 ("Ohnaka & Tsuii, 1999) up to about 90 (iMilam et al.L 
|2009|) . lEebzelter et al. (.2008., and references therein) mention 
some clustering of empirical values of '^C/'^C between 50 and 
70 for galactic C-rich stars, while their evolutionary models 
even predict values higher than 100. The solar value is about 89 
(lAnders & Grevesselll989l) . 

Line strength ratios of '^CO and '^CO transitions with iden- 
tical J provide a first order estimate of the isotope abundance 
ratio in case both lines are optically thin (Schoier & Olofsson, 
l2OO0h . In case of optically thick lines, these estimates are 
only tentative, and lower than the actual isotopic ratio. 



ISchoier & OlofssonI (l2000l) mention that flie ^^C/'^C values de- 
rived from their radiative transfer models are mostly higher than 
those estimated from intensity ratios. This effect is largest for 
their high M-objects, where the difference amounts to a factor 6. 

We note here that the line intensity ratios for 7 = 2-1 are 
consistently larger than those for 7 = 3-2 for our sample stars. 
A linear least-squares fit for the targets with line intensity ratios 
available for both 7 = 2-1 and 7 = 3-2 shows that 



/('^CO(3 - 2)) 
/(i3CO(3 -2)) 



: 4.3 + 0.48 X 



/(^^CO(2 - 1)) 
/(i3CO(2- I))' 



The standard deviation from this linear fit is 2.1 units in 
/('2C0(3 - 2))//(i3CO(3 - 2)). The slope and offset in this rela- 
tion clearly indicate that individual line intensity ratios will not 
lead to accurate values of the actual '^C/'^C abundance ratio and 
can only give rough estimates. 

We have calculated '^C/'^'C-values for our sample targets, 
correcting for line strength differences and — for the sake of 
simplicity — assuming the ta rgets were unresolved by th e beam, 
in the way this was done by ISchoier & OlofssonI (l2000h and as 
is given by Eq. [15] 



'2c/'3c 



/mb("CO(7 ^ 7 



/mb('3CO(7 -^7-1)) \vi3co(y^y-i) 



The frequency correction factor is equal to 0.87 for all values 
of 7 considered. If multiple line intensity ratios were available 
for one 7, these were averaged. The line intensity ratios and the 
resulting values for '^C/'^C are listed in Table[8] Fig. [T5] shows 
a histogram plot of the '^C/'^C-values for the AGB stars in the 
sample. It is clear from this graph that we derive very low values 
for '^C/'-'C, but as mentioned before, these estimates are lower 
limits in case of large mass-loss rate objects due to optical depth 
effects, that are not taken into account using this approach. Since 
high-M objects are more easily observed, the sample is biased 
towards the higher mass-loss rates. The values for the 12 Miras 
range from 3.6 (carbon-rich LLPeg) up to 14.4 (oxygen-rich 
RCas) and for the 3 OH/IR stars they range from 4.9 (OH 44.8- 
2.3) up to 7.2 (AFGL 5379). The spread on '^c/i^c for the other 
stars in the sample is quite large and we do not see any clear dif- 
ferences in the first-order estimates of '^C/'^C when considering 
the different types of stars in the sample. 

pVlilametal. (2009) found '^c/'^C -13 for NMLCyg and 
'^C/'^C -14 for IRC +10420, while for VYCMa they derived 
a much higher value between 25 and 46. We find similar ra- 
tios with a first order estimate for VYCMa of 14.8, and values 
of 8.7 and 8.1 for NMLCyg and IRC +10420, respectively. We 
note that Milam et al. (2009) deconvolved the line profiles for 
VYCMa and VHya into separate velocity components, since 
the different struct ures probed could have discrepant chemical 
properties (see also lZiurvs et al.L l2007. for VYCMa). They anal- 
ysed the multiple components for '^C/'^C separately, resulting 
in three different values for VY C Ma (46 , 33 an d 25) and 2 val- 
ues for VHya (71 and 72). Mil am et al.1 (l2009b claim a +30% 
uncertainty on their '^C/'-'C-values. 

5.3. Wind driving efficiency 

lOlofsson et al] (Il993h and lOudmaiier etalj (l2008h suggest fliat 
the acceleration of winds around C-rich stars is more efficient 
than that of winds around O-rich stars because of differences 
in dust absorption properties. Comparin g samples of 40 S-type , 
77 oxygen-rich and 61 carbon-rich stars. iRamstedt et aP (l2006h 
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also find indications for a higher terminal velocity in carbon-rich 
stars. 

Fig. [16] shows histograms of Vco in our sample for the 
three different chemical types of objects: oxygen-rich, carbon- 
rich and S-type stars. The mean Vco-values for respectively 24 
oxygen-rich and 13 carbon-rich AGB targets in our sample 
are 14.5 km s"' and 15.4k ms~'. Thes e values are ~4-10kms"' 
higher than those of Ramstedt et al.l (2006). We note that the 
largest V po-values for our carbo n-rich and S-type sample exceed 
those of iRamstedt et al.l (l2006l) by about 10km s '. These ex- 
treme values are however the only ones in both samples — 
V Hya in the C-rich and Gru in the S-type sample — so they 
are only weakly reflected in the samples' mean values. Our re- 
sults do not lead to the same clear conclusion that the wind 
driving efficiency is higher for carbon-rich than for oxygen-rich 
stars. They do, however, point in the same direction consider- 
ing the sampling bias we have in favour of oxygen-rich objects 
and the occurrence of higher Voo-values for carbon-rich than for 
oxygen-rich targets in our sample. 

6. Conclusions 

The data set presented in this paper covers a very large number 
of high-quality measurements of CO rotational line transitions 
7 — » 7 - 1 with J ranging from 1 up to 7. Observations of both 
'^CO and '^CO are included. It is the first time that a CO data 
set allowing one for such a detailed sampling of circumstellar 
envelopes of evolved objects is presented in the literature. 

A grid of radiative transfer models calculated with 
GASTRoNOoM covering AGB stellar and wind parameters was 
used to derive a formalism allowing one to make estimates of 
mass-loss rates. The M-estimator is valid for J up to 16, while 
until now, only estimators valid for 7 up to 4 have been presented 
in the literature. We also extended the validity range of the M- 
estimator by correcting for saturation of the CO line intensities 
and therefore accounting for optically thick winds. The intrin- 
sic spread of the estimator ensures good M-estimates when all 
input parameters are well determined. We like to point out that, 
when only one transition is used to estimate a mass-loss rate, the 
7 = 2-1 transition will give the most reliable results. 

For many of the sample AGB stars, the estimates based on 
the different rotational lines result in constant M-values. The es- 
timates for 7 targets do however hint at a variability of the mass- 
loss rate, by showing a large spread (more than an order of mag- 
nitude in M) and clear trends. Uncertainties linked to the imple- 
mented H20-cooling could produce the effect of seemingly de- 
creasing M in time, but there are two convincing cases in which 
these uncertainties cannot account for the variations seen in the 
estimates. This is the case for LLPeg and WXPsc, and we have 
shown that the M-estimates for the latter target reflect the mass- 
loss histo ry constructed throu gh detailed radiative-transfer mod- 
elling bv iDecin et alJ (l2007h . The fact that we (1) do not have 
full sampUng of all rotational 7-level^ for all sample targets and 
that (2) did not perform full radiative-transfer modelling for the 
sample, of course limit our ability to prove conclusively that a 
star has experienced a constant mass loss or mass-loss modu- 
lations. This kind of detailed modelling will be presented in a 
forthcoming paper by Lombaert et al. (in preparation). The ex- 
amples discussed in this paper show, however, that the use of the 

* The rotational transitions C0(7 — > 7 - 1) observable with ground- 
based telescopes are limited to 7 = 1-0, 2-1,3-2, 4-3, 6-5 and 
7-6. Many of the higher-frequency transitions will be observed with the 
HIFI and PACS instruments on board the Herschel Space Observatory. 



estimator allows one not only to establish good mass-loss-rate 
estimates in case of constant mass-loss targets, but also to trace 
strong variability in the mass-loss rate of AGB stars, given a suf- 
ficiently large number of CO line transitions. The fact that our 
estimates bear quite close resemblance to the mass loss derived 
through detailed radiative transfer modelling, gives confidence 
that our estimator is of good quality. 

Further observations of mostly high-excitation rotational 
lines, such as 7 = 6 - 5 and 7 = 7-6 would make it possi- 
ble to make a better distinction between constant or (possibly) 
variable mass loss. The 7 = 5-4 transition will be measured by 
HIFI for many stars in our sample. These additional measure- 
ments will provide a more complete sampling of the circumstel- 
lar envelopes and the mass-loss history of these evolved stars. 
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(a) Histogram of the M-values for all AGB stars, grouped 
per pulsational type. 
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(b) Histogram of the M-values for all AGB stars, grouped per 
chemistry type. A'c is the number of carbon-rich targets, A^o 
is the number of oxygen-rich targets, and A's is the number 
of S-type targets. 
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(c) Histogram of the M-values of all AGB stars. 

Fig. 13. Histograms containing the M-values estimated from 
Eq. |9] and given in Table Q The adopted bin size is A log(M) - 
0.5, which corresponds to about a factor three difference in mass- 
loss rate between consecutive bins. 



Fig. 14. Mean estimates for M versus period of pulsation for the 
AGB stars in the sample are plotted in black. Different pulsa- 
tional types are indicated with different plotting symbols: + for 
Miras, * for OH/IR stars, V for SRa and o for SRb type stars. The 
red * OH/IR data points were obtained from He & Chen 
priv. comm.). A linear relation fits the Miras, OH/IR stars and 
semi-regular variables in our sample with periods shorter than 
850 days. M levels off to a constant value of 2.6 x 10"^ Mq yr"' 
for longer periods. The darker shaded region represents a factor 
three spread around this fit, the paler grey region represents a 
spread of a factor 10. See Sect. l5.1.3l for further details. 
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Fig. 15. Histogram with empirical estimates of the '^C/'^C iso- 
tope abundance ratio for the AGB stars in our sample with an 
adopted bin size of two units. The targets are grouped per charac- 
teristic chemistry: carbon-rich objects (top panel), oxygen-rich 
targets (middle panel), and S-type targets (bottom panel). 
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Fig. 16. Comparative histogram plots of the expansion veloci- 
ties (Voo) of the AGB stars in the sample. The bin size was set 
to Avoo = 2.5 km s"'. The targets are grouped per characteristic 
chemistry: carbon-rich objects (top panel), oxygen-rich targets 
(middle panel), and S-type targets (bottom panel). 
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Appendix A: CO data - figures and tables 

All '^CO and '^'CO data of the sample and the soft-parabola fits 
(Sect. 12.3.1b to these data are shown in Figs. IA.11 and lA.2l The 
most important line parameters — main-beam temperature at the 
line centre and integrated intensity — and the /3 parameter from 
the soft-parabola fitting procedure as described in Sect. l2.3.Tl are 
given in Tables [ATI and IA3] References for, or remarks on, the 
data are listed in Tables|A2]and|A3]for the '^CO and '^CO data 
respectively. 



Table A.l. Overview of all data presented in Fig. lA. II The objects are ordered according to their IRAS numbers and are specified with their velocity in the frame of the Local Standard of Rest (vlsr) 
and the terminal velocity (Vo<,) of the CO envelope. For each measured transition, the main beam brightness temperature at the centre of the line (rMB.c) is given in K, together with the integrated 
intensity (/mb). given in Kkms"'. For the high-quality data, a soft parabola fit (see text) was performed and the according /J-values are listed as the third entry per transition. In case of multiple 
observations of a specific transition, the data are given in the exact same order as in Fig. lA. 11 
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1.58 




1 1 1 

1 AO 

0.69 


1 1 nn 

1 J.UU 

1 7 09 
i / .uz 

6.35 


S S 1 

J.J 1 

S 89 
J.oZ 

A 11 
O.J 1 


4 01 

1 A4 
1 .04 


41 40 

1 4 89 
14". oZ 


7 7A 
/ . /O 

O 01 

y.uj 


- 


- 




A 44 


S 1 A8 
J 1 .Oo 


1 1 SS 

1 J.J J 








13428-6232 


GLMP 363 










- 








00 


9 41 










- 


- 


- 








i J'+DZ-ZoU/ 


w riya 


41.0 


8.5 














9 A1 
Z.Oi 


"^A SO 
jO. Jw 


O 44 
U.44 


1 OA 
J.UO 


^1 .Uo 


O OS 

u.yj 








3.06 


32.06 


9 91 
Z.Zl 


1 il9 1 Q4-9SSS 


KA DOU 


2.0 


9.0 


1.10 


16.47 


U.J J 


1 1 

J. ly 


sn 1 9 

jU. IZ 


n 1A 

U.JO 


1 0"? 

j.yj 


S8 00 

jo.uy 


AO 
U.OU 








O A'X 
U.4J 


4 1 1 
1 J 


4 AO 






















9 OA 

z.uo 


in 07 
ju.y / 


\J. IKJ 


1 S/l 
l.j4 


91 10 
Zj. jU 


O 7S 
U. / J 








9 S7 
Z.J / 


10 09 
jU.UZ 


1 OS 
j.Uj 








1^1 Q/l ^11^ 


TT T iir, 

11 ljUp 


-15.0 


23.0 














4.66 


144.41 


1 A4 
1 .04 


4.93 


155.90 


1 A1 
1.01 








4.89 151.43 


1 77 
1 . / / 






















4.55 


144.38 


1.52 


3.96 


128.39 


1.40 








3.18 


82.09 


3.28 


16262-2619 


Alpha Sco 


















0.02 


0.46 






















16269+4159 


GHer 


20.4 


13.0 








0.41 


6.15 


4.29 


1.33 


23.21 


3.30 








1.18 


17.02 


2.88 








17123+1107 


V438 Oph 


9.8 


4.2 








0.20 


1.19 


1.99 


0.31 


1.60 


1.75 




















17150-3224 


Cotton Candy Nebula 


15.0 


24.0 














0.78 


18.84 


5.15 




















17411-3154 


AFGL 5379 


-21.2 


25.0 








1.59 


49.95 


2.65 


2.81 
1.76 


80.36 
50.45 


3.41 
3.48 


2.36 


64.80 


3.38 








1.83 


54.44 


3.00 


17443-2949 


PN RPZM 39 


















0.24 


5.26 






















17501-2656 


V4201 Sgr 












0.52 


1.00 




0.30 


0.11 






















18050-2213 


VXSgr 


6.5 


24.3 








0.74 


32.47 


0.30 


2.41 


94.47 


0.73 


1.18 


42.12 


0.95 














18059-3211 


Gomez Nebula 


2.7 


6.0 








0.47 


2.79 


3.86 


0.91 


4.37 


8.95 





















Continued on Next Page. . . 



Table A.I. Continued 



IRAS Alternative name vlsr \ CO 1-0 CO 2-1 CO 3-2 CO 4-3 CO 6-5 CO 7-6 







km/s 








Tmb.c 


4lB 




^MB.c 






Tmb.c 


^MB 




^MB.c 






^MB.c 


^MB 


y3 


18100-1915 


OH 11.52-0.58 


















-0.60 


-3.27 




0.02 


1.91 










-0.66 


0.07 




18257-1000 


V441 Set 


1 15.0 
















0.23 


2.94 




0.13 


4.76 










-0.06 


-7.52 




18308-0503 


AFGL 5502 


43.0 


14.7 








43.96 


253.44 40.33 


45.72 


289.44 


33.42 


52.00 


301.50 


38.30 


16.27 


94.75 34.52 








18327-0715 


OH 24.69 +0.24 


40.0 




- 


- 


- 








-0.09 


-7.40 




-0.10 


-2.95 










0.91 


-1.46 


- 


18333+0533 


NXSer 




17 8 








1.20 


26.08 


2.72 


1.29 


28.15 


3.16 


0.97 


22.02 


1.96 


0.20 


A f\C\ 


z.4y 








18348-0526 


OH 26.5+0.6 


27.4 


17.0 








-0.18 


/.U2 




1.09 


23.57 


2.19 


0.85 


18.98 


1.85 














18361-0647 


OH 25.50 -0.29 


27.5 
















1.76 


3.29 




2.31 


6.62 










0.05 


1.57 




18397+1738 


IRC +20370 


-0 5 


16 


3.65 


84.14 


1.12 


4.52 




z.U / 


2.93 


60.34 


2.21 


3.45 


66.09 


2.85 


- 






3.10 


49.77 


3.51 


18432-0149 


V1360 Aql 


40.0 
















-0.41 


6.38 




-1.23 


-2.43 










1.47 


-5.88 




18448-0545 


RSct 


57.0 


6.0 












- 


0.77 


3.85 


8.17 










- 


- 








18460-0254 


V1362 Aql 


100.0 


„ 








0.32 


lo.l / 




-0.35 
-0.30 


6.81 
-7.24 




0.02 


-5.16 










0.16 


1.21 




18476-0758 


S Set 


15.0 


8.0 


0.13 


1.27 




0.12 






0.20 


1.54 






















18488-0107 


VI 363 Aql 


75.0 








- 








2.21 


18.75 


- 


0.64 


9.98 


- 


- 






0.01 


-0.29 


- 


18498-0017 


V1365 Aql 


60.0 


„ 








0.35 


lU. /z 




0.36 
0.07 


13.07 
3.99 




0.34 


8.50 










-0.12 


1.11 




19067+0811 


VI 368 Aql 


60.0 












- 


- 


-0.62 


1.50 




-0.51 


0.44 






- 


- 


0.01 


5.61 




19110+1045 


KJK G45.07 


60 


25 




















32.42 


729.07 


6.42 














19114+0002 


VI 427 Aql 


99 


40 














1.69 


84.28 


2.47 




















19126-0708 


W Aql 


-25.0 


20.0 


0.97 


27.03 


0.93 


2.76 


OA 00 


1 1 ^ 

i.lD 


3.89 


96.71 


2.51 


4.70 


119.33 


2.32 








4.97 


134.55 


1.49 








0.99 
1.27 


26.93 
34.06 


1.71 
1.77 


1.88 


Di.jcS 


Oil 


4.75 
5.02 


118.44 
127.23 


2.51 
2.26 


4.06 
4.08 


97.28 
96.31 


2.21 
1.07 


- 






1.73 


29.21 


0.31 






„ 


„ 














2.16 


52.98 


2.49 




















19192+0922 


OH 44.8-2.3 


-70.5 


17.7 








0.38 


y.o / 


1 /IT 

1.4Z 


0.63 


15.40 


1.99 


041 


9.99 


1.50 


0.53 


y. jy 


2. J / 








19244+1115 


IRC +10420 


75.5 


42.3 


0.90 


49.94 


2.10 


1.61 


94.76 


1.44 


3.23 


178.76 


1.93 


2.94 


149.88 


2.72 


2.77 


139.62 


2.46 






















2.46 


158.99 


0.87 


3.75 


227.95 


1.73 


2.10 


118.49 


1.81 


1.63 


80.74 


2.08 








19283+1944 


AFGL 2403 


29.0 


17.0 








0.54 


10.48 


3.01 


0.24 


5.53 


2.30 


0.40 


8.83 


2.43 














19486+3247 


Chi Cyg 


10.5 

„ 


8.5 

„ 


3.20 
0.80 


35.55 
10.01 


0.38 
1.06 


6.32 
4.11 


90.86 
57.77 


0.26 
0.60 


3.49 
5.32 
5.46 


47.70 
75.11 
75.53 


0.75 
0.61 
0.62 


10.25 


125.99 


1.30 


1.45 


23.08 


0.14 








20075-6005 


XPav 


-18.5 


11.0 














1.24 


20.18 


1.14 




















20077-0625 


IRC -10529 


-18 


16.5 














1.91 


36.97 


3.27 


1.56 


33.48 


2.13 








0.49 


19.74 




20120-4433 


RZ Sgr 


-30 


9 














1.13 


12.76 


1.92 


0.80 


9.78 


1.99 








0.03 


-1.88 




20396+4757 


VCyg 


15 


15 








2.65 


43.61 


2.62 


























21419+5832 


Mu Cep 


23.0 


35.0 








0.04 


2.63 




0.23 


13.46 


0.85 


0.44 


28.31 


-0.20 














21439-0226 


EP Aqr 


-33.5 


11.5 








5.12 


30.07 75.47 


7.88 


42.15 72.53 


7.94 


42.19 68.50 


7.33 


31.35 99.32 








21554+6204 


GLMP 1048 


-18.7 


18.6 


_ 


_ 


_ 


0.53 


12.23 


2.74 


0.54 


10.22 


3.82 


0.26 


4.42 


5.65 


0.07 


1.20 


3.17 




_ 


_ 


22177+5936 


OH 104.9+2.4 


-26.0 


18.3 








0.20 


5.09 


1.46 


0.46 


10.56 


2.39 


0.11 


3.37 


0.61 


0.22 


6.50 


0.44 


-0.54 


-10.81 




22196-4612 


pil Gru 


-11.0 


30.0 














3.84 


99.14 


9.01 


3.32 


79.75 


10.77 








1.81 


39.85 


8.23 


23166+1655 


LL Peg 


-31.0 


16.0 


4.75 


101.41 


2.02 


6.45 


118.99 


3.23 


4.50 


76.16 


4.00 


4.85 


87.58 


3.43 


1.02 


22.40 


1.95 


7.79 


152.69 


2.73 










1.56 


26.97 


2.61 








2.40 


40.23 


3.53 


4.03 


72.19 


3.28 














23320+4316 


LP And 


-17.0 


14.0 


4.28 


83.92 


1.16 


7.90 


153.59 


1.72 


3.81 


67.60 


2.25 


4.14 


74.17 


2.34 


3.09 


54.90 


1.85 
















1.63 


33.70 


0.94 
































23558+5106 


RCas 


25.5 


13.5 








2.68 


56.45 


1.14 


5.20 


99.03 


2.08 


6.05 


109.97 


2.21 


5.01 


76.20 


3.42 









Table A.2. References f or t he presented '^CO data . In case of multip l e obse rvati ons of a specific tra nsiti on, the data are give n in the exact same order a s in Fig.lA. 1 | Referen ces or remarks are (1) 
'Knapp & Morris' (1985"), (2)'Mar2ulis et al.' (' 1990), (3)'Nvman et al.' (' 1992), (4) Sopka et al. ( 1989), (5) Knapp et al.' ('1982), (6) '01ofsson et al.' (1993), (7) Neri et alj ( Il998h . (8) archive data, (9) 
■Tevssier et aL C2006 ). (10) Groenewegen et al. (1998), (11) Huggins et al. (1988), (12) Wang et al. ( 1994), (13) Wannier & Sahai ( 1986), (14) Bieging & Latter ( 199^ (-) this study. 



IRAS 


Alternative name 


CO 1 


-0 


CO 2- 


-1 


CO 3 


-2 


CO 4 


-3 


CO 6 


-5 


CO 7-6 






Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 




NML Cyg 






JCMT 




JCMT 




JCMT 




JCMT 






00192-2020 


T Get 






JCMT 




JCMT 




JCMT 




JCMT 






01037+1219 


WX Psc 


BLT 


1 


JCMT 




JCMT 




JCMT 




JCMT 




JCMT 






FCRAO 


2 


OVRO 


5 


APEX 




APEX 








APEX 






OSO 


3 


JCMT 




APEX 




APEX 








APEX 






OSO 


4 




















01246-3248 


R Scl 


SEST 





SEST 





APEX 




APEX 








APEX 


01304+6211 


V669 Cas 






JCMT 




JCMT 




JCMT 




JCMT 






02168-0312 


o Get 


SEST 


3 


JCMT 




JCMT 




JCMT 




JCMT 






03507+1115 


IKTau 


IRAM 


7 


IRAM 


7 


APEX 




JCMT 


8 






APEX 










JCMT 


8 


JCMT 


8 


APEX 


















IRAM 




APEX 














04566+5606 


TX Gam 






JCMT 




JCMT 




JCMT 










05073+5248 


NV Aur 






JCMT 




JCMT 




JCMT 




JCMT 






05524+0723 


Alpha Ori 






JCMT 




JCMT 




JCMT 


















JCMT 


8 
















06176-1036 


Red Rectangle 










APEX 














07209-2540 


VY GMa 


SEST 


3 


JCMT 




JCMT 




JCMT 




JCMT 




JCMT 


07399-1435 


Galabash Nebula 










APEX 














09448+1139 


R Leo 


IRAM 


9 


IRAM 


9 


APEX 








CSO 


9 
















CSO 


9 












09452+1330 


GW Leo 


IRAM 


10 


IRAM 


10 


JCMT 


10 


CSO 


9 


CSO 


9 








NRAO 


1 1 


SEST 


6 


CSO 


12 
















SEST 























10131+3049 


RW LMi 


IRAM 


9 


IRAM 


9 


CSO 


9 


CSO 


9 


CSO 


9 




10197-5750 


GSG 08608-00509 










APEX 














10329-3918 


U Ant 


SEST 





SEST 





APEX 




APEX 








APEX 


10350-1307 


UHya 


SEST 


6 


SEST 


6 


APEX 




APEX 








APEX 


10491-2059 


V Hya 


SEST 





SEST 





APEX 




APEX 








APEX 


12427+4542 


Y GVn 


IRAM 


7 


IRAM 


7 






CSO 


9 


CoU 


9 




13269-2301 


RHya 


IRAM 


9 


IRAM 


9 


JCMT 








CoU 


9 












J L.M i 






n 

y 




















NRAO 


13 
















13428-6232 


GLMP 363 










APEX 














13462-2807 


W Hya 










APEX 




APEX 








APEX 


14219+2555 


RX Boo 


IRAM 


n 

y 


IRAM 


9 


JCMT 








JCMT 














JCMT 




CSO 


9 






CSO 


9 




15194-5115 


II Lup 










APEX 




APEX 








APEX 














APEX 




APEX 








APEX 


16262-2619 


Alpha Sco 










JCMT 














16269+4159 


GHer 






JCMT 




JCMT 








JCMT 




_ 


17123+1107 


V438 Oph 






JCMT 




JCMT 














17150-3224 


Cotton Candy Nebula 










APEX 














17411-3154 


AFGL 5379 






JCMT 




JCMT 




APEX 








APEX 














APEX 














17443-2949 


PN RPZM 39 










JCMT 














17501-2656 


V4201 Sgr 






JCMT 




JCMT 














18050-2213 


VX Sgr 






JCMT 




JCMT 




JCMT 










18059-3211 


Gomez Nebula 






JCMT 




JCMT 














18100-1915 


OH 11.52-0.58 










APEX 




APEX 








APEX 
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Table A.2. Continued 



IRAS 


Alternative name 


CO 1 


-0 


C0 2- 


1 


CO 3- 


-2 


CO 4-3 


C0 6- 


-5 


CO 7-6 






Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 


18257-1000 


V441 Set 










APEX 




APEX 






APEX 


18308-0503 


AFGL 5502 






JCMT 




JCMT 




JCMT 


JCMT 






18327-0715 


OH 24.69 +0.24 










APEX 




APEX 






APEX 


18333+0533 


NX Ser 






JCMT 




JCMT 




JCMT 


JCMT 






18348-0526 


OH 26.5+0.6 






JCMT 




JCMT 




JCMT 








18361-0647 


OH 25.50 -0.29 










APEX 




APEX 






APEX 


18397+1738 


IRC +20370 


IRAM 




IRAM 


7 


APEX 




APEX 






APEX 


18432-0149 


V1360 Aql 










APEX 




APEX 






APEX 


18448-0545 


RSct 










APEX 












18460-0254 


VI 362 Aql 






JCMT 




JCMT 




APEX 






APEX 














APEX 












18476-0758 


S Set 


SEST 


5 


SEST 


5 


APEX 












18488-0107 


VI 363 Aql 










APEX 




APEX 






APEX 


18498-0017 


VI 365 Aql 






JCMT 




JCMT 




APEX 






APEX 














APEX 












19067+0811 


VI 368 Aql 










APEX 




APEX 






APEX 


19110+1045 


KIK G45.07 














JCMT 








19114+0002 


VI 427 Aql 










APEX 












19126-0708 


WAql 


NRAO 


14 


JCMT 




APEX 




APEX 






APEX 




SEST 


3 


SEST 




APEX 




APEX 






APEX 






SEST 








APEX 




APEX 




















APEX 












19192+0922 


OH 44.8-2.3 






JCMT 




JCMT 




JCMT 


JCMT 






19244+1115 


IRC +10420 


IRAM 


7 


JCMT 




JCMT 




JCMT 


JCMT 














IRAM 


7 


JCMT 


8 


JCMT 8 


CSO 


9 




19283+1944 


AFGL 2403 






JCMT 




JCMT 




JCMT 








19486+3247 


ChiCyg 


IRAM 


7 


IRAM 


7 


CSO 


9 


JCMT 8 


CSO 


9 








NRAO 


14 


JCMT 


8 


JCMT 


8 






















JCMT 


8 










20075-6005 


XPav 










APEX 












20077-0625 


IRC -10529 










APEX 




APEX 






APEX 


20120-4433 


RZ Sgr 










APEX 




APEX 






APEX 


20396+4757 


VCyg 






JCMT 
















21419+5832 


Mu Cep 






JCMT 




JCMT 




JCMT 








21439-0226 


EP Aqr 






JCMT 




JCMT 




JCMT 


JCMT 






21554+6204 


GLMP 1048 






JCMT 




JCMT 




JCMT 


JCMT 






22177+5936 


OH 104.9+2.4 






JCMT 




JCMT 




JCMT 


JCMT 




JCMT 


22196-4612 


pil Gru 










APEX 




APEX 






APEX 


23166+1655 


LLPeg 


IRAM 


9 


IRAM 


9 


APEX 




APEX 


CSO 


9 


APEX 




OSO 


3 






CSO 


9 


CSO 9 








23320+4316 


LP And 


IRAM 


7 


IRAM 


7 


JCMT 


8 


JCMT 8 


CSO 


9 








OSO 


3 


















23558+5106 


RCas 






JCMT 




JCMT 




JCMT 


JCMT 
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Fig. lA.ll (continued) 
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Fig. lA.ll (continued) 
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Fig. lA.ll (continued) 
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Fig. lA.ll (continued) 




Fig. lA.ll (continued) 
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Fig. lA.ll (continued) 
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Fig. lA.ll (continued) 



Table A.3. Same as TablelAl] for "CO-data presented in Fig.FOl 



IRAS 



Alternative name vlsr 



CO 2-1 



CO 3-2 



'3 CO 4-3 



CO 6-5 







km/s 


^MB.c 




P 


7"mb,c 


^MB 


P 


7"mb,c ^mb P 


Tmb.c Imb P 


- 


NML Cyg 


2.0 


33.0 


0.21 


9.49 


1.52 


- 


- 


- 




_ 


00192-2020 


TCet 


23.0 


7.0 


-0.05 


0.34 


- 


0.03 


0.50 


-1 .41 






01037+1219 


WXPsc 


9.0 


19.8 


0.36 


12.93 


0.07 


0.48 


16.01 


0.48 




_ 






" 


" 


0.23 


8.09 


0.11 


0.34 


11.39 


0.21 






7J 






" 


- 


- 


- 


0.36 


11.44 


0.50 






01246-3248 


RScl 


-17.0 


17.0 


- 


- 


- 


0.11 


5.17 


-1.55 






02168-0312 


oCet 


46.4 


8.1 


1.39 


5.22 


19.70 


3.13 


13.35 


16.91 




5.17 26.54 11.62 


03507-H1115 


IKTau 


33.8 


18.5 


0.28 


10.34 


-0.48 


0.34 


13.79 


-0.68 






04566-H5606 


TXCam 


10.8 


21.2 


0.10 


3.72 


-0.17 


0.29 


10.88 


0.06 






05073+5248 


NV Aur 


3.0 


19.2 


0.22 


6.75 


0.58 


- 


- 


- 






05524+0723 


Alpha Ori 


3.5 


14.0 


0.09 


2.46 


0.03 


- 


- 


- 




_ 


07209-2540 


VY CMa 


21.0 


46.5 


0.04 


3.32 


0.57 


0.18 


12.01 


1.07 






10491-2059 


VHya 


-17.0 


30.0 


- 


- 


- 


0.14 


3.43 


8.49 






13269-2301 


RHya 


-10.0 


12.5 


0.47 


2.78 


5.26 


- 


- 


- 






14219+2555 


RXBoo 


2.0 


9.0 


0.14 


2.58 


-1.32 


0.57 


5.51 


3.04 






15194-5115 


IILup 


-15.0 


23.0 


- 


- 


- 


0.76 


31.06 


-0.06 








" 


" 


" 


- 


- 


- 


0.92 


32.48 


0.54 






17411-3154 


AFGL 5379 


-21.2 


25.0 


0.16 


6.05 


0.84 


0.05 


2.70 


- 






18050-2213 


VX Sgr 


6.5 


24.3 


- 


- 


- 


0.11 


4.29 


0.67 






18059-3211 


Gomez Nebula 


2.7 


6.0 


- 


- 


- 


0.23 


1.81 


3.25 






18308-0503 


AFGL 5502 


43.0 


14.7 


22.97 


77.97 


102.80 


24.38 


95.37 


87.94 






18397+1738 


IRC +20370 


-0.5 


16.0 


- 


- 


- 


0.28 


7.41 


0.76 






19114+0002 


V1427 Aql 


99.0 


40.0 


- 


- 


- 


0.43 


22.20 


2.12 






19126-0708 


W Aql 


-25.0 


20.0 


0.14 


5.03 


-0.26 


0.23 


8.10 


-0.28 










55 


" 


- 


- 


- 


0.25 


8.68 


0.49 








" 


" 


" 


- 


- 


- 


0.23 


7.68 


0.48 




_ 


19192+0922 


OH 44.8-2.3 


-70.5 


17.7 


0.10 


2.56 


0.13 


0.08 


2.04 


0.49 


- 


- 


19244+1115 


IRC +10420 


75.5 


42.3 


0.17 


10.57 


1.17 


0.37 


23.49 


0.93 






19486+3247 


Chi Cyg 


10.5 


8.5 








0.24 


4.65 


-0.89 






20396+4757 


VCyg 


15.0 


15.0 


0.16 


3.16 


0.20 












21439-0226 


EP Aqr 


-33.5 


11.5 


0.67 


3.18 


208.64 


1.42 


6.76 


136.35 


0.87 3.44 16.11 




22196-4612 


pil Gru 


-11.0 


30.0 








0.10 


2.81 


5.25 






23166+1655 


LLPeg 


-31.0 


16.0 








0.86 


18.36 


1.62 






23558+5106 


RCas 


25.5 


13.5 


0.11 


2.71 


0.50 


0.37 


8.01 


1.03 







b 

UJ 



o 
o 



c 

o 

o 

ft 
p- 



p. 



O 

s 

re, 
I 

a. 



Table A.4. References for the presented '^^CO data. In case of multiple observations of a specific transition, the data are given in the exact same order as in Fig. lA.2l References or remarks are (1) 
archive data, (-) this study. 



TD A C 


Alternative name 


13 r^r\ 1 

K.\J /— 1 


13 r^r\ 1 o 


13 r^c\ A 


13 r^r\ A fi: 






Telescope Reference 


Telescope Reference 


Telescope Reference 


Telescope Reference 


- 


NML Cyg 


JCMT 


- 


- 


- 


00192-2020 


TCet 


JCMT 


JCMT 


- 


- 


01037+1219 


WX Psc 


JCMT 
JCMT 
- 


JCMT 
APEX 
APEX 


- 
- 
- 


- 
- 
- 


01246-3248 


R Scl 


- 


APEX 


- 


- 


02168-0312 


o Cet 


JCMT 


JCMT 


- 


JCMT 


03507+1115 


IKTau 


JCMT 1 


APEX 


- 


- 


04566+5606 


TX Cam 


JCMT 


JCMT 






05073+5248 


NV Aur 


JCMT 




- 


- 


05524+0723 


Alpha Ori 


JCMT 


- 


- 


- 


07209-2540 


VY CMa 


JCMT 


JCMT 


- 


- 


10491-2059 


V Hya 


- 


APEX 


- 


- 


13269-2301 


R Hya 


JCMT 


- 


- 


- 


14219+2555 


RX Boo 


JCMT 


JCMT 


- 


- 


15194-5115 


II Lup 


- 


APEX 
APEX 


- 


- 


17411-3154 


AFGL 5379 


JCMT 


JCMT 


- 


- 


18050-2213 


VX Sgr 




JCMT 






18059-3211 


Gomez Nebula 


- 


JCMT 


- 


- 


18308-0503 


AFGL 5502 


JCMT 


JCMT 


- 


- 


18397+1738 


IRC +20370 


- 


APEX 


- 


- 


19114+0002 


V1427 Aql 


- 


APEX 


- 


- 


19126-0708 


W Aql 


JCMT 
- 
- 


APEX 
APEX 
APEX 


- 
- 
- 


- 
- 
- 


19192+0922 


OH 44.8-2.3 


JCMT 


JCMT 


- 


- 


19/44+1 1 15 


IKC +ltJ4zU 


JCMI 


JCMI 






19486+3247 


Chi Cyg 




JCMT 1 






20396+4757 


VCyg 


JCMT 








21439-0226 


EP Aqr 


JCMT 


JCMT 


JCMT 




22196-4612 


pil Gru 




APEX 






23166+1655 


LL Peg 




APEX 






23558+5106 


RCas 


JCMT 


JCMT 







b 

UJ 



o 
o 



c 

o 

o 

ft 
p- 



p. 



o 

s 
re, 
I 

a. 



E. De Beck et al.: CO line survey on evolved stars and M-formulae, Online Material p 18 




--^oooooo 



[>l] ' 




o lt) o in o in 
CN — ^ — o o ci 
d d d d d d 



o 








o 


sbula 


JCMT 


> 


E 

o ^ 
o ^ 


)mez Ni 


1 = 

CN 

1 

O 
O 




o 
o 


o 


J? 


lO 






w 






CO UD 

O O O 


CN O 
O O 



[>l] ■ 



[>l] ' 






o o o o o o o 



















O ^ — ^ 












o = 
^ 



[^] ' 




[>l] ' 



= ' 




oooooooo 



CN 




O 




CJ 










— o 


d d 


d d 


[>i] 






o o o o o o 



[>l] ' 




o o o o o 



o o o o o o 




o o o o 




[>l] ■ 



^ o - 



[>l] ' 




Fig. A.2. Profiles of "CO with fitted soft-parabola profiles. 
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Fig. lA.2l (continued) 
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Appendix B: Dereddening 

Since the effective temperatures of late-type stars are well- 
determined by the intrinsic V-K colour (Sect. IC.ll ). we need 
a way to derive this colour accurately. The magnitudes of both 
the V and K bands — at 0.550//m and 2.2 yum, respectively — 
are reddened due to interstellar and circumstellar contributions. 
The data were corrected for in terstellar extinction using the dust 
maps of lSchlegel et"an ([T998h . which give the interstellar red- 
dening at infinity, Eb-v,oo, for a given set of galactic latitudes 
b and longitudes I. This value was transformed into interstellar 
reddening at a given distance d using 



8000 



Eb-V — £fl-V 



1 -exp(-10 sin \b\) 

kpc 



(B.l) 



following iFeastetai] (Il990l) . The resulting valu e for Eb-v was 
combined with Ry - 3.1 and Ak - 0.112 Ay (ISchlegel et al.L 
I199&) to give interstellar extinction coefficients for the V and K 
bands. 

Circumstellar extinction was not accounted for, even though 
it is assumed to be signi ficant for high mass-loss rate LPV tar- 
gets. iKnapp et alJ (Il998h use the mass-loss rate M to determine 
the circumstellar extinction in the K band. Since the aim of the 
paper is to establish empirical relations between stellar param- 
eters and M, this m ethod can no t be a pplied to the data for the 
sake of consistencv. lKnapp et alJ (l2003i) suggest a formalism to 
calculate the circumstellar extinction in the K band, Ak^cs, as 
a function of Sn/S 2.2, the ratio of flux densities at 12yum and 
2.2/im, re spectively . This method is only valid for low optical 
depths ( Knapp et al.l l2003h . When applied to the sample, unre- 
alistic values for A/^ cs as high as several tens (LP And) or hun- 
dreds (CW Leo) of magnitudes are indeed obtained for the ob- 
jects with presumed high M values. To decide, however, whether 
or not this correction should be applied for a given target seems 
arbitrary. For these reasons we did not correct for circumstellar 
reddening. 

Appendix C: Basic stellar parameters 

In this appendix we describe the methods that we used to obtain 
the efl'ective temperature, stellar luminosity and distance towards 
objects. The results are given in Table|6]in the paper. 

C.1. Effective temperature 

Since broad-band near-infrared colours are very good tempera- 
ture criteria, being relatively independent of stellar gravities and 
abundances, they are often used to determine effective temper- 
atures. The V-K colour is a good Teff-indicator for cool, late- 
type stars, e.g . GKM gi ants and supergia nts (di Benedetto, 1993; 
iBesseU et all 1998; Be rgeat et al.Ll200lHLevesque et all 120051) . 
For hotter A t hroug h K stars the V - I colour is more suitable 
(lBesseUetal.Lll998l) . 

In all Teff versus V-K relations the Johnson-Cousins-Glass 
photometric system is used. K magnitudes for the sample stars 
that were obtained in the systems of 2MASS (K^) and SAAO 
(KsAAo) were transformed via 



K^Ks + 0.044 = 7i:sAAO - 0.0275. 



(C.l) 



following (iBessellL l2005h . Since the accuracy of the Ksaao- 
measurements (-'0.03 mag; [OHvier et al., 2001) is almost a fac- 
tor 10 better than that of ^Tg-data (-0.2 mag), preference was 
given to the SAAO-measurements when both were available. 



5000 



4000 



3000 



Bessell et al,, 1998: K and M giants, 1.5^(V-K)o^ 10 
dl Benedetto, 1993: G-K (ll-lll-IV-V) 
di Benedetto, 1993: F-G-K-M (l) (0.65V-K) 
di Benedetto, 1993: M type (IV-V), Miras, Carbon stars 
Bergeat et al., 2001: Carbon stars, 4.99^(V-K)o^ 7.0 
Bergeot et al., 2001: Carbon stars, 7.0^CV-Kjo^ 8.31 
Levesque et oL, 2005: Red Superglonts, 2.9^(V-K)q^ 8-0 



(V-K)o 

Fig. C.l. Different calibration relations linking the (dereddened) 
broad-band colour ( y-A')n and T^fe. Rela tions were taken from 
'di Benedetto* ('1993');'Be ssell et al.l dT998l) ; lBergeat et all (l2001h ; 
Levesque et al.. (.2005.) . 

The errors on the magnitudes are expected to be small relative 
to the amplitude in the K band caused by pulsations of the stars. 
The error bars on the Teff values determined here, were calcu- 
lated assuming an error on {V-K)() of 1 mag, leading to rather 
conservative temperature estimates. 

In all equations listed {V-K)o is the intrinsic V-K broad- 
band colour with (1) both the V and K magnitudes at maximum 
light, or (2) both at their mean value if maximum light was not 
available, and (3) corrected for interstellar reddening in the way 
described in Sect. |B] 

C.1.1. Long-period variables: Miras, semi-regulars and 
irregular variables 

iBessell et al.l ( Il998l) defined a T^ff versus V-K calibration for 
KM giants (luminosity class III), given by 

Teff = 9037.597 - 3101.282 x(y-ii:)o (C.2) 
+717.7044 x(y-/:)^ 
-85.83809 x(y-7i:)^ 
H-5.021194x(y-7i:)^ 
-0.1137841 x(y-/:)f,. 

This relation is valid for 1.5 mag < (V-K)() < 10 mag (see 
Fig. lC.lb . However, Bes sell et al.l ( Il998l) mention that their static 
model atmospheres cannot adequately represent the M giants 
cooler than 3000 K, since most of these are long-period variables 
(LPVs) with shock waves traversing and greatly extending their 
atmospheres. 

di Benedett^ (Il993h presented a temperature calibration in 
which Teff is derived from angular-diameter measurements based 
on Michelson interferometry. In case of variable stars, measure- 
ments at maximum light were used. The effective temperature of 
GK stars of luminosity classes II-III-IV-V is described by 

logTeff = fl-(/7-H0.0118)(y-/:)o (C.3) 

+o.om{v-K)l 

for 1.4 mag < {V-K)o < 3.7 mag, for FGKM supergiants by 

log Teff = (fl - 0.205) -(b- Om6)(V-K)o (C.4) 
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Table C.l. Effective temperatu re T^ff versus spectr al type for 
oxygen-rich targets as given bv iMarigo et alJ (12008 *). They are 
used in this paper in case of oxygen-rich LPVs which have their 
dereddened {V-K)o broad-band colour outside the valid ranges 
for the listed Tefj versus {V-K)o relations. 



Spectral Type 


MO Ml 


M2 M3 


M4 


M5 


r.ff (K) 


3850 3750 


3650 3550 


3490 


3397 


Spectral Type 


M6 M7 


MS M9 


MIO 




T,g (K) 


3297 3129 


2890 2667 


2500 





Table C.2. Coefficients a and b for the Teff calibration from V- K 
colours via Eqs.|C3]and|C4l Taken from ldi Benedettol(ll993h . 



Spectral type Luminosity class 


(V- 


-K)o 


a b 


G-K 


II-III-IV-V 


1.4 


-3.7 


3.927 0.122 


M 


III 


3.7 


- 10 


3.833 0.101 


M 


II 


3.7 


- 10 


3.859 0.108 


F-G-K-M 


I 


0.6 


- 10 


3.954 0.133 


M 


IV-V 


3.7 


- 10 


3.954 0.133 


M 


Mira 


3.7 


- 10 


3.954 0.133 


C 


Carbon star 


3.7 


- 10 


3.954 0.133 



for (V—K)o > 0.6 mag and for M-type giants of type IV, Miras 
and carbon stars by Eg. IC.4I for 3.7mag< {V-K)o <10mag 
(see Fig. IC.ll l. The constants a and b in Eqs. IC.3l and lC.4l and 
the validity ranges in terms of (V-K)o are listed in Table lC.2] 

For oxygen-rich LPVs, except for OH/IR stars, t he eff ective 
temperature was derived via Eg. IC.4l(ldi Benedettolll993h . The 
uncertainty on the temperature would then be of the order of 1 % 
according to di Benede tto ( 1993i) . 

iBergeat et al.l(l2001l) focussed on carbon-rich objects and de- 
rived a relation for Teff versus V-K given by Egs. IC.5I IC.6I 
and IC.7I The notation Teff d denotes that the relation is based 
on angular-diameter measurements, Teff is the final temperature 
adopted for their sample stars. 

logTeff.d = -0.079(rf.oo8)(V-^)o + 3.91(±o.o2) (C.5) 

for 4.99 < (V-K)o < 7 

logTeff.d = -0.061(rf.oo8)(V-^:)o + 3.79(±o.o2) (C.6) 

for7 < iV-K)o < 8.31 

logreff,d = 1.003 log Teff- 0.009 (C.7) 

These T^ff values are e stimated to be accurate to within an er- 
ror margin of +140K dBergeat et al.L 1200? ). Because this Teff- 
calibration is based on a sample of 54 carbon-rich stars, pref- 
erence was given t o this calibration, rather than to the one by 
Idi Benedetto! (1 1993 -). in case of carbon stars. The validity range 
was adopted to be 4.99 < {V-K)o < 8.31 mag, c onsistent wit h 
the minimum and maximum in Table 5 of iBergeat efaH (l2001l) . 

It is clearly visible from Fig. IC. II that the carbon stars have 
lower predicted effective temperatures than t he oxygen-ric h 
stars. This is in agreement with the results of Marigo (2002), 
showing a clear dichotomy in Teff between both chemical types 
of AGB stars. This dichotomy is caused by different molecular 
opacities in stellar atmospheres with differing C/O abundance 
ratios. 

When using these Teff estimators, one should be aware of 
the effect of circumstellar reddening on the V and K band data. 
This reddening is stronger for objects with higher valu es of M 
and depend s on the chemistry of the central star (Knapp et al.L 
Il998l l2003h . As mentioned in Appendix E the data were not 
corrected for circumstellar reddening. 

For oxygen-rich targets with {V-K)o outside the valid ranges 
for the above mentioned relations, the temperature correspond- 
ing to the sp e ctral t ype listed in Table |6] was used (Table IC.ll 
iMarigo et al.1 l2008h . An error bar of 500 K was assumed for 
these values. For the carbon-rich targets with {V-K)o outside the 
valid ranges the effective temperatures found in literature were 
used. References are given in Table |6] 



C.1.2. OH/IR stars 

In the case of OH/IR stars the photosphere is not easily probed 
since the CSE is optically thick and conseguently the infrared 
sources often do not have optical counterparts. Therefore, the 
V-K colour is no longer a meaningful diagnostic of Tss- For 
these objects an effective temperature of 2750 + 750 K was as- 
sumed. 

0.1.3. Red supergiants 

The reff-calibration of iBessell et alj ( Il998h is not valid for 
red supergiants (R SGs) because of their low logg-values. 
iLevesgue et alJ (12005) derived a Teff-scale for RSGs in the Milky 
Way given by Eg. IC.8I This relation between the broad-band 
colour V-K and the effective temperature is consistent with the 
results from their grid of marcs models and is valid for 2.9 mag < 
{V-K)o < 8.0 mag. Adopting this new Teff-calibration for RSGs 
yields a better agreement between evolutionary tracks and obser- 
vationally determine d locations of RSGs in the HR-diagram (see 
'Massev et all |2008[ F ig. 3). This T eff-scale for RSGs closely 
follows the relation of iBessell et alJ (U998) defined for KM gi- 
ants (see Fig. IC.ll ), but since Eg. IC. 81 was derived especially for 
RSGs, preference was given to this relation over Eg. IC.2I It is 
given by 

Teff = 7741.9- 1831.83(y-/:)o (C.8) 
+263.135(V-K)l - 13.1943{V-K)1. 

iLevesgue et alJ (i2005h expect that the effective temperatures for 
M-type supergiants were obtained with a precision of 50 K. 

For RSGs with (V-K)q < 2.9 or (V-K)q > 8.0, the tem- 
perature was determined using the scale of Tefr versus spectral 
type — in case of galactic targets — presented bv iLevesgue et alJ 
(120051) . 

0.1 .4. Yellow hypergiants 

The adopted Teff-cahbration of idi Benedettol (il993h for FGKM 
targets of lumi nosity class I s eem s in accordance with the results 
presented by de Jager (1998") and Oudmaii er et al] (i2008i) for the 
yellow hypergiants (YHGs) IRC H- 10420 and AFGL2343. Since 
both objects are known to exhibit large Tsff- va riations on time 
scales of a few decades dOudmaiier et al. I l2008h . this calibration 
is somewhat uncertain. 

0.1.5. Alternative methods 

iHaniff etal.i (i 19951) derive effective temperatures from angular 
diameters and bolometric fluxes obtained through fits to optical 
and infrar ed photometry. The infrared flux method (IRFM), pre- 
sented by iBlackwell & Shalhsl (il977i) . consists of determining 
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O-rlch and S-type MIros (Feast et qL, 
C-rich MIras (Feast et al., 
-rich and C-rlch Miros > F.:;a'-;t .-: nl., 
C-rich MIros (Groeneweqcr - I r I 
lOOOd^P^ 2000d (VVt.i ' ' 



996) 



-I-/- 0.39 mag rms scatiei t;;i ■/;H/ir-; sources 
0-rich and S-type SRs (a & b) 
(Yesllyoprak & Aslon, 2004 + Whltelock et al., 2000) 
SRb variables (Barthes et oL, 1999 - Grojp 3) 



i 500 



Period [days] 



(1199^: 



Yesil vaprak & Asian 
Whitelock et al 



Fig. C .2. P-Lj, relatio ns t aken f rom|Groenewegen & Whitelocld 
iFeast eta .1 (Il989l) : IWhitelock et aL (,19911) : 
2004 : iBarthes et al.l (1 19991) . and 
Dashed lines indicate extrapola- 
tions of the P-Li, relations from the literature to shorter or 
longer periods of pulsation. See text for details. 

the ratio of the bolometric flux to the flux in a selected pho- 
tometric band and comparing these results to predictions from 
model atmospheres. Since the goal of this paper is not to individ- 
ually model all sample targets, we did not pursue these methods, 
nor did we perform full fits of the spectral energy distributions 
(SEDs). The latter would require detailed radiative transfer anal- 
ysis for each star 

C.2. Luminosity 

The luminosities of Miras, OH/IR stars, semi-regular 
pulsators of types a and b and red supergiants were 
determined through empirical period-bolom etric mag- 
nitud e relations presente d in t he liter ature (Fea st et al 
19891; fOroenewegen & Whitelockl 1 1 9961: : Whiteloc k et al 



1991 



Barthes et all I1999E IWhitelock & FeastL I200C : 

Yesilvaprak & AslaJ 12004 . These relations are then used 
to calculate luminosities adopting the solar bolometric magni- 
tude Mboi,0 = 4.75 mag. 

C.2.1 . Miras and OH/IR stars 



iFeast et aP (Il989l) derived period-luminosity, P-L*, relations 
for oxygen-rich and carbon-rich Miras in the LMC with periods 
rangin g from 115 up to 420 days. Groenewegen & Whitelock 
(Il996h revised the f-L* relation for carbon-rich Miras 
and extended the per iod range up to periods of 520 days. 
IWhitelock et al.l (Il991h developed a similar relation for OH/IR 
targets with periods of pulsation exceeding 1000 days. These re- 
lations - and their extrapolations to shorter or longer pulsation 

periods - are shown in Fig. IC.2I 

The relation derived bv iFeast et al.l (Il989l) for oxygen-rich 
Miras was assumed to be vaUd also for galactic oxygen-rich 
Miras and was adjusted with a distance modulus of 18.5 mag. 
It is given by 



IWTiitelock et al.l (Il99lh mention a rather large scatter of 
0.39 mag on their derived f-L* relation for OH/IR variables 
with periods ranging up to 2000 days. The large range in pre- 
dicted luminosities associated with this scatter has been made 
visible in the figure. Note th at the extrapolated relation for 
oxygen-rich Miras derived by iFeast et al] (1 19891) matches the 
higher luminosity edge derived for OH/IR targets very well. 
IWhitelock et al. ( 1991) also point out that there are suggestions 
that the luminosities of OH/IR sources with P > 1000 days 
probably follow this extrapolate d P-Lj, rela t ion th at was estab- 
lished for oxygen-rich Miras bv iFeast et al.l (Il989l) . Because of 
these arguments and the strong evolutionary connection between 
oxygen-rich Mh' as and OH/IR star s (see Sect. 12.1b we adopt the 
P-Li, relation of lFeast et aLld 19891) for oxygen-rich Miras, given 
in Eq. IC.9I for the determination of luminosities for all oxygen- 
rich Miras and OH/IR targets, regardless of their period of pul- 
sation. Relation [C.9| was also used in determining the luminosity 
of S-type Miras in the sample. 

Luminosities of carbon-rich Miras were determined via 



Mbol,C-MIRA = -2.59 log P + 2.02(±o.26) 



(C.IO) 



following iGroenewegen & White loc5 (Il996l) . This relation pre- 
dicts lum inosities substantia lly higher than does the relation de- 
rived by iFeast et al] (|1989|). Preference w as given to the for- 
mer since the sample of iGroenewegen & Whitelock (1996) is 
almost three times as large as the one used by F east et al.l (119891) 
(54 versus 20 objects). For longer-period carbon-rich Miras, this 
P-Li, relation was extrapolated, similar to what was done in the 
oxygen-rich case. 

C.2.2. Semi-regular pulsators: SRa & SRb 

Many semi-regular type a (SRa) objects diff'er from Miras only 
by their smaller amplitudes in the V band (AV < 2.5 mag). 
Objects classified as semi-regular b (SRb) variables have rather 
poorly defined periodicity, show a superposition of multiple pul- 
sation periods, or have alternating intervals of periodic and slow 
irregular changes. 

Semi-regulars and Miras are linked in terms of evolution, 
since S Rs of types a and b are thought to be the progenitors o f 
Miras (IWhitelock & FeastL 120001: lYe^ilvaprak & AslanL |2004 . 
Therefore some type of f -L* relation, analogous to those for 
Mi ras is probably justified for th ese LPVs. 

lYe^ilvaprak & Aslaiil (|2004 presented period-magnitude re- 
lations for all magnitudes in the Johnson system and the IRAS 
system for SRa and SRb stars. To determine the bolometric mag- 
nitude, and hence the luminosity, only the J and K period - mag- 
nitude relations were used. 

J = -1.03(±o.56) log/' -3.47(^0.98) (C.ll) 

K = -1.38(±o.58)logP- 3.99(^1.00) (C.12) 

^ Yesilvaprak & AsianI mention a scatter of 0.63 mag and 0.85 mag 
for the J and K relations, respectively. The bolometric correction 
on the K magnitude was calculated via 



/Wbol,0-MIRA,OH/IR - -3.00±o.24 log P "H 2.85±o.f 



(C.9) 



BCk = 10.86- 38.10(7 -/:)o 

-h64.16(7 - K)l - 50.72(7 - K)l 

-Hl9.48(7 - K)^o - 2.94(7 - K)l 
taken from lWhitelock et al.l (l2000l) . where 



(C.13) 



with P in days. 



Mboi = K + BCk. 



(C.14) 



E. De Beck et al.: CO line survey on evolved stars and M-formulae, Online Material p 23 



Relations IC. 11 1 and IC . 1 2 1 were however only adopted in es- 
timating the luminosities for SRa variables, and not for the 
SRb variables in the sample. Barthes et al. (1999) presented 
P -{J - K) relations for long-period varia bles, including a lar ge 
number of SRb type stars. The sample of Barthes et al.l ( 1 999h is 
divided in four groups, and the SRb type stars are largely split 
up in two of those — Groups 2 and 5, corresponding to younger 
and older kinematic properties respectively. Their Group 2 is far 
less luminous than their Group 3 and preference is given to the 
relations for their Group 3 - Eqs. IC.lSl and IC.16I - since these 
estimates are in best accordance with luminosities derived in the 
literature for the SRb objects in the sample. Bolometric correc- 
tions were again derived via Eq. IC.13I 

/TsRb = -2.11(i0.i3)(logi^-1.75)-6.51(±,,o5) (C.15) 
(J-K)o = 0. 17(^0.04) log P + 0.85(^0.08) (C.16) 

The uncertainties on calculated luminosities only take into 
account the uncertainties on the coefficients in Eqs. lC.9l through 
IC.16I The error bars on the pulsation periods are assumed to 
be small compared to the errors on the coefficients in the above 
relations. 



C.2.3. Red supergiants 

The luminosity of the RSGs in the sample was determined 
as a function of period of pulsation and effective temperature. 
Eq. lC.17| gives the relation between the K magnitude and the pul- 
sation period f or RSGs in the milky way galaxy as presented by 
lDambisl(ll993h . This relation was completed wit h the bolometric 
correc tion on the K band magnitudes given by iLevesque et al.l 
(l2005h (Eq. ICTSl l. This method leads to luminosities ranging 
from ~10"' Lq up to ~5 10^ Lq. It is given by 



OH/IR stars, since these provide a very accurate distance es- 
timate ("van Langeveld e et al.L Il990h . (2) Hipparcos-data were 
used when available with a relative error lower than 50%. (3) 
values in the literature derived from SED fitting were used and 
assumed to have uncertainties of 50 %, unless stated otherwise 
in the original papers. 

C.3.1. Alternative method 

An alternative method to derive distances is to compare ob- 
served magnitudes with absolute magnitudes predicted through 
formalisms as presented in Sect. lC.2l The relations that could be 
used for distance determination from magnitudes for Miras are 



^O-MIRA = -3.47(±o.I9)logP + 0.98(±0.45) 

^^c-MiRA = -3.561ogP+ 1.14 



(C.19) 
(C.20) 



Since the K magnitude is less sensitive to metallicity than 
the bolometric magnitude mboi, it is a better distance indicator 
dWhitelock & Feasn, l2000i and references therein). The strong 
variabiUty in the K band of most sample stars, however, gives 
rise to large uncertainties on absolute magnitudes. An additional 
constraint to the application of this method to determine dis- 
tances is that the circumstellar material can redden the K band 
magnitudes significantly. 



K = -10.10ao.08)- 3.44(logP-2.7) 



BCk = 5.574- 0.7589 X 



1000 K 



(C.17) 
(C.18) 



C.2.4. Yellow hypergiants 

No general period-luminosity relations can be established for 
these targets because of their very strong variability. The latter is 
reflected in the mul tiple motions from blue to red and vice versa 
in the HR diagram dde Jaged [1991 lOudmaiier et al.L 12008 '). For 
both AFGL 2343 and IRC +10 420, die two YH Gs in the sample, 
luminosities were taken from Ide Jaged (Il998h . yielding values 
for log L*/L0 of 5.30 and 5.80, respectively. 



C.2.5. Alternative methods 

More recentlv. lBm-thes et alJ (Il999l) and lKnapp etaP (l2003h de- 
rived P - K relations for oxygen-rich LPVs, including Miras 
and semi-regulars of types a and b. Using these relations in 
determining the luminosities of Miras would introduce extra 
uncertainties, since bolometric corrections for K band magni- 
tudes would have to be derived. This is the main reason why 
the above mentioned P-Mboi relations of' Feast et al.l (119891) and 
[Groenewegen & Whitelock (1996.) were used for Miras. 



C.3. Distance 



Distances towards the sample stars have been determined in var- 
ious ways. (1) OH-maser phase lag distances were used for most 



